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1.0 ABSTRACT 


This method uses a panel method to calculate incompressible flow about 
arbitrary three-dimensional inlets with or without centerbodies for four 
fundamental flow conditions: unit onset flows parallel to each of the coord- 

inate axes plus static operation. The computing time is scarcely longer than 
for a single solution. A linear superposition of these solutions quite rigor- 
ously gives incompressible flow about the inlet for any angle of attack, 
angle of yaw, and mass flow rate. Compressibility is accounted for by apply- 
ing a well-proven correction to the incompressible flow. Since the computing 
times for the combination and the compressibility correction are small, flows 
at a large number of inlet operating conditions are obtained rather cheaply. 
Geometric input is aided by an automatic generating program. A number of 
graphical output features are provided to aid the user, including surface 
streamline tracing and automatic generation of curves of constant pressure, 

Mach number, and flow inclination at selected inlet cross sections. This 
report describes the inlet method, including the use of the program and presents 
illustrative results. 
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3.0 GENERAL DESCRIPTION OF THE METHOD 

The computer program that has been constructed for analyzing flow in and 
about the forward portions of nonaxi symmetric engine inlets applies to nonround 
inlets, inlets with curved centerlines, and scoop inlets. The range of operating 
conditions to which the method is applicable is essentially unrestricted except 
that the flow at large distances must be subsonic with respect to the inlet. 
Generally speaking, the present method may be considered a generalization of 
that of Stockman (Ref. 1), which is restricted to axisymmetric geometries. 

No attempt has been made to develop new physical or mathematical theories. 
Instead, well -proven computational techniques have been assembled in a new way. 
The main criterion used in assembling this method has been numerical efficiency 
with special emphasis on the case where solutions for a fairly large number of 
distinct operating conditions for a given geometry are desired. A second major 
consideration has been to assure maximum utility for users of the procedure. 

This has been accomplished by providing a geometry generating program to facili- 
tate input and certain graphical output options to display the calculated flow 
quantities in a manner shown by experience to be useful. 

3.1 The Panel Method 

The basic calculational tool is a panel method for calculating inviscid 
incompressible lifting flow about arbitrary configurations (Ref. 2). In 
recent years the use of panel methods has grown to the point where at least the 
simpler, older versions are now standard aerodynamic design tools. Moreover, 
the literature on panel methods has increased in a similar degree. Accordingly, 
details of the panel method are not given here, but only some aspects that are 
particularly important for the present application are discussed. 

The body about which flow is to be calculated is input to the computer by 
means of the coordinates of a set of points lying on the body surface. These 
points are used to form a large number of quadrilateral "panels," which form the 
discrete representation of the body. An example is shown in Figure 4. On each 
panel there is a surface source distribution whose strength is initially unknown. 
The calculational procedure forms a large full matrix that represents the 
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aerodynamic influences of the panel sources on each other. A matrix equation is 
then solved, which expresses the condition that the total flow about the body 
must have zero normal velocity at one point of each panel. That is, the source 
strengths on the panels are adjusted to cancel the normal components of what- 
ever onset flow or flow environment is incident to the body. Each right side of 
the linear equations corresponds to a particular onset flow, and yields a com- 
plete set of panel source strengths and thus a complete flow about the body. The 
aerodynamic influence matrix depends only on the geometry of the body. If, as 
in the present program, a direct matrix solution is used, solutions for several 
flows about the body may be obtained in essentially the same computing time as a 
single solution. 

Because panels are placed on the actual surface of the body, rather than on 
some interior surface, the method described above is applicable to arbitrary con?- 
figurations with no restrictions to slender bodies or small flow perturbations. 
Moreover, the method is numerically exact in the sense that any degree of accuracy 
may be obtained by using a sufficiently large number of panels. The key calcula- 
tional efficiency of a panel method, as opposed to, say, a finite-difference 
method, is that the former need not solve for the entire flow field, but can obtain 
the flow only on the surface if desired. Flow at off-body locations can then be 
generated only where it is of interest, e.g., at certain interior cross sections. 

3.2 The Technique of Superposition of Fundamental Solutions 

The basic panel-method calculational efficiency as described above is aug- 
mented in the inlet application by additional efficiencies arising from use of 
the technique of superposition of fundamental solutions to obtain flows at vari- 
ous inlet operating conditions. For a simple closed body in incompressible 
flow the operating condition is defined by the magnitude and direction of the 
flow at large distances. Thus there are three parameters: angles of attack and 

yaw, and velocity magnitude. An inlet, however, contains complicated internal 
machinery that controls the amount of fluid that enters the inlet. It is not posr 
sible or even desirable for the panel method to analyze this machinery in detail. 
Instead, its effect is lumped into a single parameter, mass flow through the inlet. 
Thus, in incompressible flow the operating condition of an inlet is defined by 
four parameters. In compressible flow there is a fifth parameter defining compres- 
sibility, e.g. sound speed or total temperature. 



To analyze the flow about the inlet, the forward portion of the inlet, 
which is the region where the flow is of interest, is artificially extended by 
means of a long afterbody with constant inner and outer cross sections, as shown 
in Figure 1(a). The afterbody is open at the downstream end. A forward location 
is selected by the user and the interior cross section there is designated the 
control station where mass flow is evaluated. Often the propeller plane or the 
compressor face is chosen as the control station. Four fundamental flow solu- 
tions are calculated. The first three are illustrated in Figure 1. Incompres- 
sible flow about the inlet with its extension is calculated by the panel method 


(a) 



Inlet geometry with afterbody. 



V 


(b) Zero angle of attack and yaw 




(c) Pure angle of attack and pure yaw 


Figure 1. Three fundamental flow solutions for three-dimensional inlets. 
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for three inclinations of the uniform onset flow at large distances: (1)0° angle 

of attack and yaw, (2) 90° angle of attack with 0° angle of yaw, and (3) 90° angle 
of yaw with 0° angle of attack. Mass flow for these three fundamental flows are evalu- 
ated at the control station. In general, these cannot be predicted in advance, but 
whatever values are computed simply must be accepted. Usually the mass flows for 
flows (2) and (3) are near zero. The above three fundamental flows may be linearly 
combined to give the flow about the inlet at any angle of attack and any angle of 
yaw, but the mass flow will always be nearly equal to that of flow (1). To con- 
trol mass flow a fourth fundamental flow is required. 

The fourth fundamental flow is the incompressible flow about the extended 
inlet in static operation, i.e., zero flow at large distance but a finite mass 
flow at the control station. This is obtained as illustrated in Figure 2 which 
also shows the N-lines of input points and the lifting strips utilized by the 
panel method (Ref. 2). A known distribution of vorticity is placed on the 
panelled inlet surface in such a way that closed vortex lines run around the 
inlet cross sections — both interior and exterior. (An axisymmetric inlet has 
ring vorticity oriented exactly on the circular cross sections.) The vorticity 
strength is constant in the direction parallel to the long afterbody, i.e., in 
the axial or stream direction. The variation of vorticity strength circumfer- 
.entially "around" the inlet, i.e., from one lifting strip to another, may be 
specified by the user if he desires, but the strength is normally taken constant 
which is definitely the standard option. Possibly better results would be 

LIFTING STRIP 



Figure 2. Use of surface vorticity to obtain a solution for a three-dimensional 
inlet in static operation. 
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obtained for an unusual inlet geometry using a "circumferentially" varying 
vorticity strength, but this is up to the user to determine. In any case, the 
normal velocity on the inlet surface due to the vorticity distribution is not 
zero, and source distributions on the panels are also required. The vorticity 
field provides one more onset flow to the inlet and one more right side of the 
normal -velocity matrix equations. The resulting solution is a fourth flow 
about the inlet. It can be shown that for a sufficiently long afterbody this 
flow represents the inlet in static operation. There is a constant uniform 
flow interior to the inlet at a sufficient distance from the entrance. On the 
exterior surface the flow velocity falls to zero with increasing distance from 
the entrance, as indeed it does in all directions in the fluid. 

Thus the same mechanism, i.e. vorticity, that is used to produce lift in 
reference 2 is used here to regulate mass flow. In the lifting application, the 
vorticity strengths on the lifting strips are determined by the Kutta condition. 

In the inlet application the distribution of vorticity is prescribed and the 
total strength is adjusted to give the desired mass flow. 

The four fundamental flow solutions may be linearly combined to give the 
incompressible flow about the inlet at any operating condition. Specifically, 
the four combination constants are determined so as to yield any prescribed 
values of: (a) angle of attack, (b) angle of yaw, (c) freestream velocity, and 

(d) mass flow at the control station. For incompressible potential flow the 
validity of this superposition is not approximate but is exact. The computational 
effort required to perform the combination for a particular set of operating con- 
ditions is a small fraction of that required for the original panel -method cal- 
culation. Thus solutions for large numbers of operating conditions may be 
obtained quite cheaply. In the inlet program provision has been made for perma- 
nent storage of the four fundamental flow solutions, so that additional solutions 
may be obtained easily at future times. Not only may additional operating con- 
ditions be specified, but the location of the control station may be changed, 
and additional interior cross sections may be designated where it is desired to 
know the flow field. 
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Since the computing time for the panel method is much less than that 
required for a finite-difference solution even for a single operating condition, 
it is estimated that if several operating conditions are of interest, the 
present program is faster than a finite-difference method by two orders of mag- 
nitude. This numerical efficiency is realized because the main flow calcula- 
tion of the present program is incompressible. Thus the practicality of this 
program depends, in an essential way, on the availability of an accurate and 
general compressibility correction that may be applied a posteriori. That is, 
it is crucial that a means be available for obtaining compressible flow about an 
inlet from the incompressible flow about the same inlet, as opposed to, say, a 
stretched version of the inlet. 

3.3 The Compressibility Correction 

The compressibility correction used in the present program is that of 
Leiblein and Stockman (Ref. 3). The mathematical details are contained in the 
references and will not be repeated here. However, a general description of the 
procedure appears to warrant inclusion. 

The Leiblein-Stockman compressibility correction is a correlation based on 
empirical observation. The correlation was deduced from a comparison of exact 
solutions for the compressible and incompressible flows in a turbine stator pas- 
sage. Strictly, the correction applies to internal flows, but it has been 
extended to external flows as well. 

Physically, the correction consists of one-dimensional compressible flow 
in the predominant flow direction. Specifically, the compressible total-to- 
static density ratio is obtained from one-dimensional flow considerations. This 
requires only the flow passage area, the inlet mass flow, and the total temper- 
ature. This one-dimensional density ratio is then modified by the ratio of local 
incompressible velocity to the average incompressible velocity at that particular 
axial station. The local incompressible velocity is obtained at each point by 
combining the four fundamental flow solutions. The average incompressible 
velocity may be obtained either from the incompressible solution or from one- 
dimensional incompressible flow. For points in the interior of the inlet, the 
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key quantity in this correlation is the flow passage area, which is automatically 
calculated by the program from the inlet geometry. At points on the exterior of 
the inlet an alternate procedure is employed. 

Of course the key issue determining the validity of using a compressibility 
correction of this sort is how well it predicts compressible flow as obtained 
experimentally. The Leiblein-Stockman correction is exceptionally well verified 
in this regard (Refs. 1, 3, 4). Indeed it appears to have attained the status 
of a routine procedure. While all previous work has been done with round inlets, 
excellent comparisons with experiment have been obtained at large angles of 
attack — at least 75° — where the flow is essentially three-dimensional in nature. 
Rather than repeat the comparisons contained in the references, a new comparison 
is presented here. Figure 3 compares calculated and experimental static-to-total 


a = 75°, M jh = 0.603 



X NONDIMENSIONAL DISTANCE FROM HIGHLIGHT 

Figure 3. Comparison of calculated and experimental static-to-total pressure 

ratio distributions on the windward side of an inlet operating condi- 
tion at which the flow is supercritical. 
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pressure ratios on the surface of an inlet at 75° angle of attack. The agree- 
ment is excellent even at the pressure peak where the local Mach number is 1.32. 
This illustrates the validity of the compressibility correction for shock-free 
flows with significant supersonic regions. 

3.4 Geometric Specification of the Inlet 

For good accuracy the surface of the inlet must be specified by a rather 
large number of points - larger than is needed for purely external flows. To 
minimize the effort required of the user of the inlet program, a geometry pack- 
age has been incorporated into it. This geometry package has a wide range of 
capabilities and is applicable to a much larger class of problems than inlets 
(Refs. 6, 7). In the present application the geometry package allows the user 
to input a relatively small number of points defining the inlet and centerbody. 
Moreover, extreme care in distributing the points need not be taken. The 
routine enriches the point number and redistributes the points according to one 
of a number of algorithms. Graphical output of the points allows the user to 
approve the distribution before proceeding with the relatively expensive flow 
calculation. Alternatively, the user may elect to enrich the point number and 
continue to complete the flow calculation in a single computer run. An example 
of the use of the geometry package is shown in Figure 4, which presents the input 
point distribution and the enriched distribution for a typical case. 

3 . 5 Symmetry 

It is anticipated that most inlets will have a plane of symmetry. Accord- 
ingly* provision has been made to take advantage of this condition to reduce 
both computing time and the amount of required input data. One half the body is 
defined by points, and the reflected half is accounted for internally. There 
are two aerodynamic influence matrices — one matrix for flows where the values 
of the source densities on symmetrically located panels are equal and another 
matrix for flows where these values are equal and of opposite sign. Of the four 
fundamental flow solutions, all use the first matrix except that corresponding to 
pure yaw, which uses the second. Here pure yaw is defined as a freestream per- 
pendicular to the inlet's symmetry plane. 
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INPUT DISTRIBUTION 


ENRICHED DISTRIBUTION 


Figure 4. Use of the geometry package to obtain input data for the inlet program. 

/ 

Even though only half the inlet is input, probably most users prefer 
graphical output for the complete inlet. Therefore, this has been made the 
standard option, both for the inlet geometry (Fig. 4) and the calculated flow 
results. Results for the other half of the inlet are generated internally, 
including cases with yaw where the combined numerical values are not symmetric. 

3.6 Graphical Output Options 

The main tabular output of the inlet program consists of local flow proper- 
ties, e.g., velocity, pressure, Mach number, etc., at one point of each of the 
panels on the inlet's surface and at points on sections across the interior of 
the inlet. To aid the user in evaluating this rather formidable amount of data, 
two graphical output options have been provided. These are the surface stream- 
line option and the cross-section option. 

The user may specify a set of points on the surface of the inlet and/or its 
centerbody through which streamlines are to be constructed. The program calcu- 
lates these streamlines — in both upstream and downstream directions if desired. 

It then graphically displays the streamline trajectories against the panel 
representation of the inlet. Viewing angles may be specified by the user, and 
a number of angles may be input simultaneously. Either orthographic (parallel) 
projection or perspective projection may be selected. A second type of graph 
is also generated by this option, namely conventional two-variable Cartesian 
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plots of pressure, Mach number, and flow angles versus distance along each 
streamline. 

An interior cross section of the inlet is defined by a point and a vector. 
The program defines a plane through the point whose normal is the given vector, 
calculates the region of this plane bounded by its intersection with the inter- 
ior of the inlet and the centerbody, and automatically distributes field points 
throughout this region. For any inlet operating condition, flow is calculated 
at the field points, and an interpolation procedure defines curves of constant 
value of various flow quantities: pressure, Mach number, flow angle. A user 

may specify several cross sections if he desires. The particular routine used 
for constructing contours of equal value is quite general. In particular a 
contour of some particular value may consist of more than one disjoint portion. 
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4.0 GENERAL REQUIREMENTS OF THE GEOMETRIC INPUT 


The present program has been built around that of references 2 and 5, 
which is designed for lifting wings and wing-fuselages. The geometric input 
requirements reflect this ancestry. Moreover, to prevent the logic of the 
Flow Passage Area Calculation from growing too complicated, certain restrictions 
on the complexity of the inlet calculation have been assumed. 

The inlet program is designed for configurations of the type illustrated 
in Figure 4: a single isolated cowl or inlet either with or without a single 

disjoint centerbody. Additional surfaces such as slats, vanes, or additional 
inlets (say an inlet within an inlet) will cause the logic of the Flow Passage 
Area Calculation and, probably also the Cross Section Calculation, to fail. 

In the logical scheme of the lifting panel method the inlet (cowl) is a 
lifting body on which surface vorticity is placed, and the centerbody is a non- 
lifting body which has no vorticity. Both bodies are specified to the computer 
by means of coordinates of a set of points that define their surfaces. These 
points are organized by means of so-called N-lines, as indicated in Figure 2 
and described in reference 2. Basically, an N-line is a curve in the body surface 
connecting successive input points. Thus points are input N-line by N-line. In 
the inlet program the N-lines are not cross sections at a given axial station but 
rather are curves at a given circumferential location that traverse the bodies in 
the axial or streamwise direction. 

Specifically, the first point of each N-line lies at the most aft (down- 
stream) location of the extended afterbody (Figures 2 and 4). Moreover this 
point is on the interior surface of the inlet. Successive input points on the 
N-line proceed forward (upstream) along the interior of the inlet, go around the 
lip, and then proceed aft (downstream) along the exterior surface to the end of 
the afterbody. The last input on the N-line is a so-called wake point. It is 
assumed that at the aft end of the afterbody the inlet centerline is straight. 

The direction from the second-to-last point on an N-line to the last (wake) 
point should be parallel to the local centerline direction. On the centerbody 
the first point on each N-line is at the nose and successive points proceed aft 
(downstream) to the end of the afterbody. Thus, inlet N-lines both begin and 
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end at the end of the afterbody, while centerbody N-lines begin at the front and 
proceed to the rear. The order of input of N-lines on both the inlet and the 
centerbody is clockwise as viewed by an observer at upstream infinity. 

Every N-line on the inlet must have the same number of points. Every N-line 
on the centerbody also must have the same number of points, but this number may 
be different from that of the inlet N-lines. It is required, however, that both 
the inlet and the centerbody have the same number of N-lines. 

If one plane of symmetry is utilized, this must be the y = 0 plane, and 
points must be input in the plane to obtain a closed surface. Thus the first 
and last N-lines of both the inlet and the centerbody must lie in the y = 0 
plane. 

In all cases the Geometry Package of references 6 and 7, which has been 
incorporated into the inlet program, is available to aid in preparing geometric 
input data. 
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5.0 OVERALL OPERATION OF THE PROGRAM 


There are basically two main modes of operation of the inlet program. The 
first begins with the geometric input data and produces the Fundamental Solutions 
of Sections 3.2 and 6.0. This is the time-consuming portion of the calculation. 
Its logical structure is illustrated in Figure 5. Once the fundamental solu- 
tions have been generated, they can be saved indefinitely to produce flow solu- 
tions at any future time. Thus the expensive part of the calculation need only 
be done once for each inlet configuration. The second mode of operation begins 
with values of the physical quantities that define an inlet operating condition 
and produces a complete flow solution for that condition. It is illustrated in 
Figure 6, particularly 6c. The first two parts of Figure 6 illustrate the Flow 
Passage Area and Cross Section Calculations that are performed at the same time. 
It is possible to run the program all the way through both modes in a single 
computer run, but this probably will not be the standard mode of operation. In 
any event several sets of operating conditions may be input simultaneously. 


GRAPHICAL 

OUTPUT 


GRAPHICAL 

OUTPUT 



TABULAR 

OUTPUT 


Figure 5. Logic of generating the fundamental flow solutions. 
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(a) Flow Passage Program 



(b) Cross-Section and Control Station Program 



(c) Combination Program 

Figure 6. Logic of calculating flow at a specified inlet operating condition 
from the fundamental flow solutions. 
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6.0 PROCEDURE FOR NONSYMMETRIC CASES 


6.1 Basic Potential -Flow Program 

This program is very similar to the ordinary lifting panel method of ref- 
erence 2. It differs: (1) in replacing the Kutta condition by a mass flow 

condition, (2) in the procedure for cases with symmetry, and (3) in the computa- 
tion and storage of individual flow solutions. 

Solutions are obtained for three uniform onset flows, one parallel to each 
coordinate axis. Thus the flow components are: (1, 0, 0), (0, 0, 1) and 

(0, 1, 0). A solution is also obtained corresponding to a uniform vorticity on 
each lifting strip. If there is no symmetry, all solutions use the same matrices 
V.. and A.., representing velocity at the ith control point due to a unit 

I J I J 

source density on the jth panel (see reference 2 for detailed definitions). 


Source density solutions are obtained for the three uniform and L non- 
uniform vorticity onset flows, where L is the number of lifting strips. These 
last are denoted a^', j = 1, 2, ..., N, k = 1, 2, ...» L. 

The L sets of source densities for the vorticity solutions are combined 
into a single set by assigning the relative vorticity strengths B^, 
k = 1, 2, ..., L. There are two options: 

fkl 

1. Constant option: All the vorticity strengths B v ' are simply set 

equal to unity. 

(kl 

2. Input option: A set of L numbers B' ' are input. 


In either case the sets of a are combined by the formula: 


L 



k=l 


j " 1» 2, ..., N 


(1) 


where the subscript j denotes values on the jth panel and N is the total 
number of panels. This set is denoted the static solution. There are 

thus four sets of source densities and four onset flow velocities: one for 

each of the three uniform flows and one for the static solution. The onset flow 
velocity at the ith control point for the static case is given by 
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( 2 ) 


L 

y( s ) = V* B (k) V (k) 

10 w. 10 

where the V-j 0 ' are the onset flows for the individual strip vorticity solutions. 

These source densities and onset flow velocities are used with the influ- 
ence coefficient matrices in the usual way to obtain on-body and off-body vel- 
ocities. The new feature is that this is done for all four individual flows. 

These velocity distributions, and source densities for the individual flows form 
the fundamental solution data set. After these solutions have been stored, the 

matrices V . . and A., are destroyed, 

■ J 1 0 

6.2 Cross Sections and Control Station 


The cross section and control station program furnishes sets of off-body 
points to the program. Say there are C cross sections specified, each with P 
points. Moreover, each cross section has its own normal vector. Velocities at 
each of the CP off-body points are calculated for each of the four onset flows. 
For each point this involves calculating a new row of the V^. matrix including 
the vorticity onset flows, and then performing summations of the form 



(3) 


using the source densities and onset flows from the fundamental 

solution data set. The complete set of: off -body points, flow velocities, 
boundary points (where the cross-section intersects the on-body N-lines), areas, 
normal vectors and mass flows at each cross section comprise the cross section 
and control station data set. The areas and mass flows are obtained by numerical 
integration. One particular cross section is specified as the control station, 
and data for it are used in the combination program. 


6.3 Combination Program (Excluding Mach Number) 

To obtain an incompressible combined solution, the basic input consists of: 
angle of attack a, angle of yaw 3, freestream velocity and mass flow 
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Q at the control station. The combined velocity (or source density) formula 
is 


V. = VJcosa cose?! 100 ) + sinav! 001 ) 


+ cosa sin&vj°^J + Bvj s ' 


( 4 ) 


The V- 's on the right are the individual fundamental solutions. The single 
constant B is evaluated from the mass flow Q at the control station by 

Bq( s ) = q _ v^tcosa cosBQ^ 100 ^ + sinaQ^ 001 ^ + cosa sineQ^ 010 ^] (5) 

where the Q's are the individual mass flows from the control station data set. 


After an incompressible combined velocity has been calculated as above, a 
compressibility correction is applied if a nonzero Mach number has been speci- 
fied. The flow passage program calculates a flow passage area for every on-body 
and cross-section point. For each point this area is used together with 
the incompressible velocity in a compressibility correction program to yield 
compressible velocity and pressure. A normal output is accomplished for each 
Mach number. These calculations are discussed in subsequent sections. 
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7.0 PROCEDURE FOR SYMMETRIC CASES 


The symmetry procedure for inlets differs in two important respects from 
the symmetry procedures for an ordinary lifting case: the use of two flow 
matrices and the generation of twice the usual amount of output. This section 
refers to the nonsymmetric description (Section 6.0) for all portions of the 
calculation that are similar. 

7.1 Input and Geometric Output 

As is done in the ordinary lifting panel method of reference 2, only the 
half-body for which y >_ 0 is input, and the tabular output of control points, 
etc., i.e., the basic geometric quantities, is for the half-body. However, the 
graphical output (if any) is for the complete body, so before entering the 
graphics program the other half of the body is generated by reflection. Each 
reflected panel is obtained from a basic (input) panel by changing the signs 
of five geometric quantities (reference 2). Each input section has its corres- 
ponding reflected section. Thus there is twice as much geometry entering the 
graphics program as was input. 

7.2 Basic Potential -Flow Program 

-V 

The major change here is that two complete V.. and A^. are calculated 
and stored: one calculated by the procedure for a plus symmetry plane and one 
for a minus symmetry plane. These correspond, respectively, to flows that are sym- 
metric and antisymmetric about the symmetry plane. These are calculated together. 
After a V.j has been calculated for both a panel and its reflection, the two 
influences are added to get the influence appropriate to a plus symmetry plane 
and subtracted to obtain the influence appropriate to a minus plane. Both are 
saved in different places, so that finally the plus and minus V^. and A. . 
matrices are stored separately. All matrices are of order N, where N is the 
number of input panels (half the body). 

The vorticity onset flows are calculated by the procedure appropriate to a 
plus symmetry plane only. They are summed in the manner of Section 6.1 to 
obtain a static onset flow. 
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The fundamental flow solutions, source density and velocity, are obtained 
in the usual way using the plus matrix for all flows except the uniform onset 
flow (0, 1, 0) which corresponds to pure yaw and using the minus matrix for 
that one flow. In particular, the plus matrix is used for the two uniform onset 
flows (1, 0, 0) and (0, 0, 1) and for the static onset flow. Each flow solution 
is only for the half-body that was input, i.e. N on-body points and off-body 
points having y ^ 0. 


7.3 Cross-Sections and Control Stations 

Off-body points, normal vectors and boundary points are calculated as in 

the nonsymmetric case but only for y > 0. When calculating velocity influences 
*> 

from the V.. matrices both plus and minus influences must be calculated. The 

' J 

minus influences (matrix row) are used with the solution for the flow (0, 1, 0) 
and the plus influences for all other fundamental solutions. The cross-section 
and control station data set is generated in the usual way, except the computed 
areas are doubled, and mass flows are doubled except for the (0, 1, 0) 
solution where the mass flow is set equal to zero. 

7.4 Combination Program (Excluding Mach Number) 

The combination constant B is determined as in the nonsymmetric case, 
equation (5). The combined incompressible velocity is obtained from the funda- 
mental flow solutions by equation (4) for the control points of the input (y >_ 0) 
panels. Another complete set of output is produced for the control points of the 
reflected panels. The geometry of the reflected panel is obtained by changing cer- 
tain signs. In particular the reflected control point coordinates are x Q , -y Q , 
z and the normal vector is r\ . -n u » n 7 » where the unsigned quantities are 
for the input panel. 

The velocity components at the reflected control point are obtained from 
equation (4) using the same fundamental flow components but with certain sign 
changes. These changes are 
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Solution 

All Other 


for (010) 

Solutions 

Sign(s) 

x-component 

y-component 

Changed 

z-component 



This holds for off-body points as well. 

The flow passage program must be altered to reflect symmetry and to compute 
flow passage area for both basic (input) and reflected panels. The flow 
passage area corresponding to the control point of a reflected panel is set 
equal to the flow passage area corresponding to the associated basic panel. 
Compressible velocities and pressures are then computed for both basic and 
reflected panels in the same way as for a nonsymmetric case. 

Thus, there is twice as much output as for a nonsymmetric case having 
the same number of input points. 
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8.0 FLOW PASSAGE AREA CALCULATION 

The compressibility correction of reference 3 requires a flow passage area 
to be associated with each point where compressible velocities are to be com- 
puted. At points interior to an axi symmetric inlet (either on or off the 
surface), the proper choice of such an area is straightforward. It is the area 
of the circular cross section of the interior of the inlet at the axial location 
of the point in question. The area may be constructed geometrically by passing 
a plane through the point in question normal to the axis of symmetry of the 
inlet and then finding the intersection of this plane with the interior surface 
of the inlet. If a centerbody is present, its intersection with the plane is 
also computed. The flow passage area is then the area contained within the 
closed curve defining the inlet intersection minus the area contained with the 
closed curve defining the centerbody intersection. 

Determination of the flow passage area for interior points of a three- 
dimensional inlet is also straightforward and indeed can be accomplished by 
essentially the above procedure provided: (1) the inlet (and centerbody) has a 

straight centerline and (2) the inlet highlight lies in a plane perpendicular 
to the centerline. In such cases the centerline plays the role of the symmetry 
axis in the construction of the previous paragraph. However, for general three- 
dimensional inlets, e.g., inlets with curved centerlines or scoop inlets, the 
proper definition of the flow passage area is not obvious. In constructing the 
present program it was decided that the program would automatically calculate 
flow passage areas for all interior points and that this calculation would be 
entirely geometrical. To accomplish this it was necessary to assume that the 
inlet highlight lies approximately in a plane, although the orientation of this 
plane is arbitrary. The centerline of the inlet is defined by the user who 
inputs a set of centerline points. 

At points exterior to the inlet (either on or off the surface) the concept 
of flow passage area becomes less significant and its definition somewhat 
arbitrary. Nevertheless, flow passage areas have been defined for points on 
the surface in a manner that has yielded reasonable compressible velocities in 
a limited number of preliminary cases. The program automatically calculates 
these areas. Selection of flow passage areas for exterior off-body points is 
currently an unsolved problem. 
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The remaining subsections of Section 8.0 present the details of the flow 
passage area calculation. 


8.1 Organization of the Input 

The restrictions on the body geometry and the order of input are set forth 
in Section 4.0. The result is a configuration as shown in Figure 7. Specifi- 
cally, there is one inlet (lifting body) and either one or no centerbody 
(nonlifting body). The N-lines on the inlet are input from lower surface 
"trailing edge" to leading edge to upper surface "trailing edge." In this 
context the lower surface is inside the inlet, and the upper surface is outside. 
Thus, for example, the top of the inlet - both inside and outside - is a single 
N-line, and the bottom — both inside and outside — is a different N-line, as 
shown in Figure 7. All N-lines on the centerbody (if any) are assumed to start 
at a certain "initial point" and to proceed downstream. Thus, in particular, 
the top and bottom of the centerbody are different N-lines. 


The centerline of the inlet is input as a set of points and normal 
vectors beginning with a certain initial point and continuing downstream 
(Figure 7). While normally the input centerline corresponds to the physical 
centerline, it need not. The input centerline is a device for calculating flow 
passage areas. In particular, the initial centerline point and vector define the 
inlet lip or highlight. 


INITIAL POINT 
OF CENTERBODY 

INITIAL POINT 



TOP N-LINE- INLET 


OTHER CENTERLINE POINTS WITH 
NORMAL VECTORS 


TOP N-LINE- CENTERBODY 


BOTTOM N-LINE-CENTERBODY 


OF CENTERLINE 


C5 


BOTTOM N-LINE- INLET 


Figure 7. Input to the program including the centerline points used to calcu 
Tate flow passage areas. 
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8.2 Intersection Program 


The basic calculational unit is a subroutine that takes a segmented line 
(N-line) and finds the point of intersection with a plane: the coordinates of 

the intersection and the number of the segment it is on. The plane is defined 
by a point and a normal vector. The routine also has the capability of comput- 
ing the area bounded by a set of intersection points obtained for several 
N-lines. 


8.3 Initial Centerline Point 

The initial centerline point and its vector are key input quantities. 
Together they define a plane whose intersection curve with the inlet defines 
inside and outside. More specifically, this plane should intersect the panels 
on each side (interior and exterior) of the inlet highlight. This is easy to 
arrange if the inlet has a highlight that lies in a plane all around the inlet. 
The plane of the initial centerline point then is parallel and slightly down- 
stream. If the inlet highlight is not plane, the plane of the initial center- 
line point should be as near the highlight as possible and should intersect 
the inlet all the way around. More specifically the plane defined by the initial 
centerline point and its normal vector must have exactly two intersection points 
on each N-line of the inlet, otherwise the calculation ceases with an error 
message. For the centerbody there is either one intersection point on each 
N-line or there are no intersection points on any N-line. The Intersection Pro- 
gram calculates the plane's intersection points with all N-lines. 

8.4 Interior and Exterior Points of the Inlet 

Points on all N-lines are classified as either interior or exterior 
and a designation must be added to the data describing each control point. On 
the inlet for each N-line the one of the two intersection points that is closer 
to the initial centerline point is designated an interior N-line point. Then, 
on the strip of panels adjacent (following) the N-line, the control point of 
the panel adjacent to the segment containing the closer intersection point is 
designated the "last interior control point." (If the k-th segment of the 
N-line contains the intersection point, the control point of the k-th element 
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of the strip is the desired one.) It and all preceeding control points of the 
strip are labelled interior. All succeeding control points of the strip are 
labelled exterior. 

If there are no intersection points with the N-lines of the center- 
body, then all control points of the centerbody are labelled interior. If there 
is an intersection point on one segment of each N-line of the centerbody, the 
control point of the panel adjacent to that segment is denoted the "first 
interior control point." It and all succeeding control points of the strip are 
labelled interior. All preceeding control points are labelled exterior. 


8.5 Initial Flow Passage Area 


For the inlet, the interior intersection points (interior N-line points) are 
collected into one list and the exterior intersection points into another list. 
(These last are the intersection points on each N-line further from the center- 
line point.) Using the intersection program, areas are computed for both these 
sets of points: the interior area A^ and the exterior A g . If there are 

intersection points of the plane with the centerbody, the area enclosed by these 
points is called A c . Define 


A 


FP 





( 6 ) 


If there are no intersection points with the centerbody then 


A 


FP 




(7) 


Next arc lengths are computed along N-lines from the initial point to the 
intersection points: the interior and the exterior points on each N-line of 

the inlet and the point (if any) on each N-line of the centerbody. For a segment 
containing an intersection, the arc length up to the first point of that segment 
is given by the method of reference 2. The additional arc length to be added 
to get the arc length of the intersection point is found by computing the usual 
straightline distance between the intersection and the first point of the segment. 
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The area A| p is associated with the arc lengths corresponding to the 
exterior intersection for all N-lines of the inlet. The area A pp is associ- 
ated with the arc lengths corresponding to the interior intersection for all 
N-lines of the 'inlet and also for all N-lines of the centerbody, if any. 

8.6 Successive Flow Passage Areas 

The procedure of Sections 8.3, 8.4, 8.5 above is repeated for all of the 
imput centerline points. A plane is passed through each centerline point having 
its normal vector identical with the centerline tangent vector. Interior and 
exterior intersection points are found on each N-line of the inlet and a single 
intersection point on each N-line of the center body (if applicable). The 
areas A pp and A| p are calculated from (6) or (7) and associated with the 
arc lengths corresponding to the intersections (interior and exterior) on the 
inlet N-lines and the centerbody N-lines as described in Section 8.5. Thus, 
tables of area versus arc length are constructed (Section 8.8). 

8.7 Possibilities for the Centerbody 



Figure 8. Possibilities of the centerbody. 
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There are four possibilities for the location of the centerbody, some of 
which require slightly different logic. 

a. No Centerbody : This case is identified by' input flag. Equation (7) 

is used for all points of the inlet, and, of course, no area table is constructed 
for the centerbody. 

b. Centerbody Abreast of Inlet Highlight : This case is identified by the 

fact that the intersections of the centerbody N-lines with the initial centerline 
plane occur on the first segments of those N-lines. Equations (6) are used for 
all intersection points, and a table of area versus arc length is prepared for 
each N-line of the centerbody using the areas of (6). Tables for the centerbody 
have the same length as those for the inlet N-lines except for the one additional 
exterior intersection on each of the latter (Section 8.8). 

c. Centerbody Downstream of Inlet Highlight : This case is identified by the 

fact that there are no intersection points of the initial centerline plane with 
any of the N-lines of the centerbody, but there is a centerbody. This situation 
may exist for the first several centerline points. Eventually, a centerline 
point is reached that is the first whose plane intersects the centerbody. Call 
this the centerbody point (not necessarily the same as the initial centerbody 
point where the N-lines start). Equations (7) are used for the intersections 
corresponding to the centerline points preceeding the centerbody point. Equa- 
tions (6) are used for the intersections corresponding to the centerbody point 

and all succeeding points. It is only for the latter cases that there are 
intersections on the centerbody. Tables of area versus arc length are prepared 
for all N-lines on the inlet and the centerbody, but the centerbody tables are 
shorter, because their entries begin with the centerbody point, while those of 
the inlet begin with the initial centerline point. 

d. Centerbody Upstream of Inlet Highlight : This case is characterized by the 
fact that there are intersections of the N-lines with the initial centerline 
plane, and these are not on the first segments of the centerbody N-lines. Just 

as in (b), equation (6) is used to compute areas for all interior intersections 
on both inlet and centerbody N-lines, but now there are additional expressions 
to assign values of flow passage area at points of the centerbody upstream of 
the inlet lip (Section 8.10). 
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8.8 Arc Length Tables. Interpolation for Flow Passage Areas at 

Interior Control Points 

When the various tasks of Sections 8.3 or 8.7 have been accomplished, flow 
passage areas App have been associated with the arc lengths corresponding to 
interior intersection points on both inlet and centerbody N-lines. For each 
N-line a table of area App versus arc length can be constructed. 

For each inlet N-line, the first entry in the table is the area A| p 
associated with the exterior intersection for the initial centerline point. 

The second entry is the area App associated with the interior intersection 
for the initial centerline point. Successive entries are the areas App for 
the interior intersections for the successive centerline points, i.e.. 

Associated 

Entry Arc Length Area 

1 S-j Aj£ p (ext. Initial centerline point) 

2 S£ App (int. initial centerline point) 

3 S 3 App (int. 2nd centerline point) 

4 S 4 App (int. 3rd centerline point) 


Last S L ast App (int. last centerline point) 

Notice this table is in decreasing magnitude of arc length. 

For each centerbody N-line the first entry is the area A pp associated with 
the first intersection point. This is obtained from the initial centerline 
point in cases b and d of Section 8.7, but from some other centerline point in 
case c. Successive entries are the areas Ap p for successive centerline points. 
This table is in increasing magnitude of arc length. 

Interpolation in arc length is performed in the above tables to obtain flow 
passage areas at the midpoints of the segments of each N-line. 
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Once the above has been accomplished, the flow passage areas for the 
control points must be obtained. This is obtained by averaging the areas for 
the midpoints on the two segments bounding the panel on . which the midpoint 
is located. Specifically, the area to be associated with the control point of 
the k-th panel of a strip is obtained by averaging the areas associated with 
the midpoints of the k-th segments of the two successive N-lines bordering the 
strip. 

The above procedure is used to obtain flow passage areas A pp at the 
interior control points of the inlet and the centerbody body. The exterior 
control points require special handling. 

8.9 Flow Passage Areas at Exterior Control Points on the Inlet 

The area A| p defined in Section 8.5 is not really a flow passage area 
but is merely a geometric area. It is used to obtain flow passage areas at exter- 
ior control points. As mentioned above, the definition of such an area is some- 
what arbitrary, but the one adopted has given encouraging results. It is import- 
ant to emphasize that even if the results at exterior points should turn out 
less accurate than desired, it would not affect the accuracy of interior points, 
for which the compressibility correction has been thoroughly verified in a large 
number of cases. 

Initially the A| p are handled in a manner similar to the internal flow 
passage areas A pp to obtain a value of Aj* p at each exterior control point. 
Additionally, however, two special values of Ap p are needed: an initial 
value A q near the inlet highlight and a final value A L far downstream on 
the afterbody. 

Key quantities for this computation come from the intersection with the 
N-lines of the inlet with the plane for the initial centerline point. For each 
N-line there are the coordinates x^ , y i , z.. of the intersection point (interior) 
and the flow passage area A pp . All four of these quantities are averaged 
between successive N-lines (a two-point average to produce an initial point 
(x 0 , y Q , z ) and an initial area A Q for the strip of panels lying between the 
two N-lines in question). Further, the coordinates of the last intersection 
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on each N-line (from the last centerline point) are averaged between successive 
N-lines to produce coordinates x L , y L » z L of the last point on the strip 
lying between the two N-lines, and a similar average of areas AjS p gives the 
final area A^. 

Also needed are the exterior geometric areas Aj£ p for the exterior 
control points. These are obtained by interpolation in a manner similar to 
that described above for interior control points. For each N-line a table 
of area versus arc length is constructed for the intersection points. The 
first entry in the table is the area A pp corresponding to the interior inter- 
section for the initial centerline point. The second entry is the area Ajjf p 
associated with the exterior intersection for the initial centerline point. 
Successive entries are the areas A£ p for the exterior intersections for 
successive centerline points, i.e.. 


Entry 

1 

2 

3 

4 


Arc Length 
S 1 


Area 

A F p (int. initial centerline point) 
A£ p (ext. initial centerline point) 
App (ext. 2nd centerline point) 

Aj* p (ext. 3rd centerline point) 


Last 


LAST 


A£ p (ext. last centerline point) 


The table is in increasing magnitude of arc length. 


Interpolation in arc length is performed in the above table to obtain 
areas A| p at the midpoints of the segments of each N-line, whose arc length 
is calculated in the manner described in Section 8.8. These midpoint areas are 
averaged between successive N-lines to obtain areas A| p at the exterior control 
points. 
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The geometric area A£p associated with each exterior control point, 
together with the control point coordinates x, y, z, is used to obtain the 
following linear dimensions needed to compute the exterior flow passage area. 



Finally, the flow passage area associated with the control point for use 
with the compressibility correction is 

A p p = [1 + D]itR[R + x - x 0 ] (9) 

Equation (9) gives a flow passage area that varies from a value near the 
inlet highlight that is approximately equal to the highlight disk area to an 
area far downstream that equals the surface area of a sphere centered at the 
highlight minus the final cross section of the inlet. This area distribution 
is approximately that through which the flow passes in the case of static 
operation. 

8.10 Flow Passage Areas at Exterior Control Points of the Centerbody 

The first step is to compute the cross-sectional area of the centerbody at 
selected "axial" locations. To do this, a set of centerline points upstream of 
the inlet highlight are input. These are used in the same way as centerline 
points as far as computing intersections and areas except that: (1) these 
points are input and used only in case d of Section 8.7, and (2) intersections 
and areas are found only on the N-lines of the centerbody. 

The resulting areas A c are interpolated in arc length to find areas at 
the control points in the same manner as in Section 8.8. Finally, the flow 
passage area to be saved with each exterior control point on the nonlifting 
centerbody is 

a fp = a 1 -a c (10) 
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where A i is the area computed in Section 8.5 for the initial interior inter- 
section. Comparing equation (6) (which applies here rather than (7)) and the 
definitions in Section 8.9 shows that 

A i = A o + A c (initial centerline point) (11) 

The A c in (11) is that for the initial centerline point. The A in (10) is 
that for the control point in question. Also, A Q is defined in Section 8.9. 

8.11 Flow Passage Areas for Off-Body Points 

There are two types of off-body points: ordinary input off-body points 

and off-body points generated by the cross-section and control station program. 

The latter are always interior points, and the flow passage area is calculated 
by that program (Section 9.0). It is anticipated that all interior off-body 
points will be so generated. Thus, input off-body points will be exterior points. 
It is not clear at this time how the flow passage areas for such points should 
be defined, and this decision has been postponed. Thus input off-body points 
should not be used in compressible cases. 
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9.0 CROSS-SECTION AND CONTROL STATION PROGRAM 

In inlet applications it is not only the flow on the body that is of inter- 
est but also the flow at field points off the body, usually in the interior of 
the inlet. Most commonly, distributions of various flow quantities across some 
interior cross sections of the inlet are desired. Of course, off-body points 
could simply be input to the program, but this is tedious for the user. Instead, 
the program generates the required off-body points at each cross section. The 
user specifies the location and orientation of each cross section where flow 
calculation is required, and the program automatically distributes off-body 
points across the cross section. It is this last calculation that distinguishes 
a cross-section calculation from the flow passage area calculation of Section 8.0. 
One cross section is selected as the control station for the combination program 
and thus at least one cross section is required if a combined flow is to be 
computed. The remaining subsections of Section 9.0 give the details of the 
cross-section procedure. 


9.1 Input 

Each cross section is defined by a cross-section point and normal vector 
(six numbers) which are used in a manner very similar to the way a centerline 
point and its normal vector are used in the flow passage program (Figure 7). 
However, more calculations are required for a cross section. One additional 
input integer designates which cross section is to be used as the control sta 
tion in the combination program. Another input integer specifies the number 
of "radial" segments to be used at each cross section. This latter number 
is denoted r. 


9.2 N-Line Intersections and Flow Passage Area 

The intersection program is used to find the intersection points of the 
N-lines with the plane through the cross-section point whose normal vector is 
the input cross-section normal vector. On the N-lines of the inlet only the 
interior intersection points are needed. The exterior intersection points are 
discarded. If there are intersections with the N-lines of the centerbody, these 
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points are needed. The areas A. and A (if any) are computed as in 

I V 

Section 8.5, and the flow passage area is computed from (6) or (7) as appropriate. 

9.3 Cross-Section N-lines 

Two situations must be distinguished: the case when there are intersections 
with the N-lines of the centerbody and the case when there are no such inter- 
sections. 


9.3.1 No centerbody. — Cross-section N-lines are to be constructed, each 
of which starts at the cross-section point and ends at an intersection point on 
an inlet N-line. Clearly one cross-section N-line is constructed for each inlet 
N-line. For each cross-section N-line the first point is the cross-section point, 
the last is the intersection point on the inlet N-line, and the interior points 
are obtained by interpolation (Figure 9). 


Let x . y„„, z„ p be the coordinates of the input cross-section point and 
(kl (kl cs (k) cs cs 

x: , y> , z) ' be the coordinates of the intersection point on the k-th 
inlet N-line, k = 1, 2, ..., T, where T is the total number of inlet N-lines. 
The points on the k-th cross-section N-line have the coordinates 



INTERIOR INTERSECTION 
POINTS ON THE INLET 
N-LINES 



INTERPOLATED POINTS 


Figure 9. Generation of cross-section points and panels in the absence of a 
centerbody. 
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where r is the input integer mentioned in Section 9.1. Each cross-section 
N-line has (r + 1) points. 

9.3.2 With centerbody . — Cross-section N-lines are to be constructed, each 
of which starts at an intersection point on a centerbody N-line and ends at an 
intersection point on an inlet N-line. For each cross-section N-line the first 
point is the centerbody intersection point, the last is the inlet intersection 
point, and the interior points are obtained by interpolation (Figure 10). 

It is assumed that there are the same number of inlet and centerbody N-lines, 
say T of each; clearly there are T cross-section N-lines. 



Figure 10. Generation of cross-section points and panels in the presence of a 
centerbody. 
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(k} Ik) fkl 

Let > z q be the coordinates of the intersection on the k-th 

centerbody N-line and, as before, let xj k ^, y! k \ zj k ^ be the coordinates of 
the interior intersection point on the k-th inlet N-line where k=l,2, ...,T. 
The points of the k-th cross-section N-line have the coordinates 
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As before, each N-line has r + 1 points. 


9.4 Cross-Section Panels and Control Points 


Cross-section panels are formed from the points of the cross-section N-lines 
in a manner logically similar to (but using different formulas than) that used to form 
surface panels from the on-body N-lines. In particular, a panel is formed from 
four points, which consist of two consecutive points on an N-line and the cor- 
responding two points on the adjacent N-line. Symbolically, the m-th panel of 
the k-th cross-section strip is formed from the points: 
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x m-lJ 

M*r>- 

-y (k) l 

J m-lJ 

3 + 


II 

r x (w) 

LVi 

-x (k) l 

x m J 


f * 

~>f 

3 + 



(14) 


(15) 


are formed. These are in the cross-section plane. 
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The unit normal vector is the same for all panels on all strips, 
input cross-section normal vector divided by its own length, say 


" = V + v + n z R 


where 


'x N * 


"y = /- "z-JT- N = 


and where N , N , N are the input components. The unit tangent vectors are 
x y z 

the same for all panels on all strips. Define 


*2 = n x j = -n z i + Oj + n x k , 


-j- . ? 2 

T| 5 T 2 x n = ” n x n y"* + n x + V J “ n y n z k 


These vectors are divided by their own lengths to obtain unit vectors 

t = t = -1Z T = — T = JjZ + T^ + T 2 

T lx T-j ’ r ly T ] * ‘lz T ] * 'l Vlx ‘ly 1 Iz 


and the same for t 2x » t 2y , t 2z « These are used as a transformation matrix 
(just one matrix for all panels) 


a ll = ^x* 


a 12 t ly* a 13 = *lz 


a 21 t 2x > a 22 hy* a 23 l Zz 


a 31 V a 32 n y* a 33 " n z 


The area of a panel is given by 


1 -*■ -> 
A = I - R 1 x R 2 


R, 

lx 

R iy 

R lz 

TV 

ro 

X 

R 2y 

R 2z 



n 

n„ 

n 



X 

y 

z 


The control -point of the panel, which is treated as an off-body point in 
the flow calculation, has coordinates 
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Also desired are the cross-section coordinates, 5, n» of the control point. 
These are obtained by a transformation using the above transformation matrix and 
the cross-section point as origin. From formulas (18) - (20) it can be verified 
that if the cross section is normal to the x-axis, then the £ axis is parallel 
to the y-axis and the n-axis to the z-axis. In all cases the n-axis is paral- 
lel to the y = 0 plane. 


The quantities needed for each panel are 


V y o’ V n ’ A 


These are computed and stored for all panels on the cross-section — a total of 
(T — 1 )r panels. 


Quantities needed for the cross section as a whole are 


x cs* y cs* z cs 

a ll 

a 12 

a 13 

App 



a 21 

a 22 

a 23 


(24) 


a 31 

a 32 

a 33 




Also needed for each cross section are the reference and cross-section 
coordinates of the boundary points. These are just the intersection points on 
the inlet N- lines. That is, 





k = 1 , 2 T 


(25) 


f k 1 ( k 1 

together with the coordinates ^ obtained by transforming these into 

cross-section coordinates. The same is done for the centerbody intersection 
points (if any) ; 
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( 26 ) 



and their transformed cross-section coordinates 



n 


(k) 


9.5 Flow Calculation 


The data generated above are used with: the basic potential -flow program, 
the geometric quantities data set, and the source densities from the fundamental 
solution data set. Each cross-section control point is treated as an off-body 
point, and the flow there is computed. 

For each cross-section control point a row of V.j matrix is computed by 
the usual formulas. Also calculated are the vorticity velocities at the off-body 
points due to each of the individual vorticity flows (one for each lifting strip). 
These last are combined with the B's in the same manner as for the on-body 
points to yield the static onset flow velocity in the form of equation (2). 


Fundamental flow velocities are computed in a manner analogous to equation 
(3) for the cross-section control points, i.e., 


*( 100 ) 

v i 

y(010) 


■S 

J-l 


■Ev"'5 


j=l 




(27) 


V. (static) 


j=l 


V. .o 
U J 


(s) + J(S) 


It is understood that in symmetry cases a different matrix V.j is used for 
the 010 flow. 
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The above velocities are computed for each cross-section control point to 
form the cross-section data set. Also in this data set are the vector F and 
scalar F volume fluxes — one set for each fundamental solution. These are 
given by 

f- £ ^ A. (28) 

cross 

section 

and 

F = F • n (29) 

where V. is identified with each of the four velocities above. is the ith 

panel area. The scalar flux at the control station is needed in the combination 
program. In symmetry cases the F's calculated from thecbove equation are 
doubled and then the F corresponding to the 010 onset flow is set equal to 
zero. 
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10.0 COMBINATION PROGRAM INPUT AND PRELIMINARY CALCULATIONS 


10.1 Incompressible Option 

If the flow is incompressible, this option is selected and only the fol- 
lowing quantities are input: 

freestream velocity 

V c average velocity at the control station 

V ref reference velocity used in computing pressure coefficient 
a angle of attack 

3 angle of yaw 

k integer denoting that the k-th input cross section is to be used 

as control station 

The volume flux Q at the control station is obtained from 

Q = V c • App(k) (30) 

where App(k) is the area of the control station. This is used in equation (5) 
to evaluate the combination constant B. Then velocity is calculated from equa- 
tion (4). 


10.2 Freestream Conditions 

For compressible flow the freestream conditions are defined by inputting 
angle of attack a, angle of yaw 3, and three additional quantities: 

either velocity or Mach number M^ 

either total pressure P t or static pressure P g 

either total temperature T t or static temperature T $ 

Then the preliminary calculations are as follows. 
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10.2.1 M„ input. — 


(a) If P t is input, P g is given by 


p .- p - 


tO + j "ij 


If P $ is input, P t is given by 


v -. K *.) 3 ' 1 


If neither is input, the default is 


P t = 2116.23 


and P s is as above. 


(b) If T. is input, T s is given by 

T s ■ T t(’ + H)' 1 
If T s is input, T t is given by 

T t* T .(’ + H) 

If neither is input, the default is 

T t = 518.67 


and T s is as above. 

In either case stagnation and freestream sound speeds a t and a $ are 
calculated from 


a t - 49/T^, 


a s = 49/r T 


and V from 


». * a t M ~( 1 + f *i)' 
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( 39 ) 


10.2.2 Vo Input . - 

(a) If T t is input, a t is given by 

a t = 49^ 

M is then calculated from 

oo 



The remainder of the calculation proceeds as in 10.2.1 above. 

(b) If T s is input, a $ is given by 

a s = 49 /fj (41) 

and M by 

co v 

M = V /a„ (42) 

oo oo £ 

The remainder of the calculation proceeds as in 10.2.1 above. 


10.2.3 Additional Freestream Quantities . - Built into the program are 
constants 

g = 32.174 
R = 1715.63 

The following quantities are calculated: 

P t 

Total density: P t = 

P s 

Static density: p $ = 

Dynamic pressure: = 0. 

Pressure ratio: P /P* 
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Density ratio: 


p s /p t 


Temperature/sea-level ratio: e = 


Pressure/sea-level ratio: 


6 = 


518.67 

p t 

2116.23 


Critical speed: a* = a^/Zl .2 


Maximum velocity: V = /5 a. 

max t 


Equivalent Incompressible Freestream velocity: V' = V / — | 

40 ”\ p t/ 

Equivalent Incompressible Critical velocity: a; = 0.6339a* 

10.2.4 Summary . - Three freestream conditions are input: \l a or 

Pt or P s , T t or T s (or default values). Calculated and saved are 

V co’ M co> P f P S* T t* T S* a t’ a S 

p f v v y p t> p s /p t- 01 6 

«*• V max’ v :> «J 
Nineteen quantities altogether. 


10.3 Control Station Conditions 


Input consists of one of the following three quantities: 
w inlet mass flow rate 
V average velocity 

L 

M c average Mach number 

The remaining two must be calculated plus some additional quantities. 



10.3.1 V Input. - p is given by 

V W 


w by 


and PI by 
c 


p c = p t 


i 

5 ' a t 


2n 


2.5 


App(k) 

W = 9P C V C -TO- 


V. 

c “ a. 


M_ = -£ 




1 - 1/2 


(45) 


(46) 


(47) 


10.3.2 M Input. - V is given by 

v L 

v c = a t H c( 1 + F M c)" 1/2 < 48 > 

Then p and w are obtained as in 10.3.1 above. 

10.3.3 ft Input . — Here V must be calculated iteratively by solving the 
equation 


V 


c 


ft 

g(A Fp (k)/144)p t [l - l/5(V c /a t ) 2 ] 2 - 5 


(49) 


starting with V c = 0 

Once V is known, PI and p„ are obtained as in 10.3.1 above, 
c c c 

10.3.4 Additional Control -Station Quantities . — These are calculated as 
follows: 


Pressure ratio: 
Density ratio: 


<p c /p t > - 


, _i/Vf 

5 \h 


l3.5 


p c /p t 


(50) 
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Dynamic pressure: q 


= 0 - 7P t(^ 


FT / M c 


Velocity ratio: V oo /V c 

Jq 

Corrected mass flow: ft CQr = ft — 


( 50 ) 


Equivalent Incompressible Average velocity: = V c I 


10.3.5 Summary . - One quantity, ft, V , or M c is input. Quantities 
saved are 

*• v c’ H C> V (P c /P t>’ <Oc /p t>> V "/ c l' *cor’ r c 
A total of ten quantities. 


10.4 Combination Calculation 

The basic combined incompressible velocity at each on-body, off-body and 
cross-section point is obtained from equations (5) and (4), where 

Q = V'A Fp (k) (51) 

and V' is used instead of V in these equations. Once the vector F. is 

CO 00 I 

calculated from (4), the scalar magnitude is calculated by the usual square 
root of sum of squares. 

For the present discussion the i in V i denotes incompressible rather 
than the usual meaning of velocity at i-th control point. However, it is 
obvious that all calculations must be done for every control point individually, 
as well as for every off-body point and cross-section point. 

10.5 Compressibility Correction 

There is a flow passage area Ap p associated with every on-body control 
point, off-body point, and cross-section point (Section 8.0). This is now 
used in the compressibility correction. 
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First the average incompressible velocity is calculated from 


t v FP' 


If V. > a*, V. is set equal to a* and the average density ratio e is 


e = fi- 


= 0.6339 


If this occurs the point is labelled "choked" on the output. 

If V. < a;, it is used as it stands to compute (p7p t ) by an iterative 
procedure. The iterative equation is 

L ,/MM 2 - 5 


with initial value e = 1. 


e = 1 - 


5 \ a t/7J 


Finally, the compressible velocity magnitude V is calculated from 


v 


This is used to calculate all output quantities (Section 11.0). First V is 

tested. If V > V m . set V = V mav and use M = ~ and P. / P . = 0 instead 
— max max s' t 

of values calculated below, and denote this point as "vacuum." If V < V _ 

max 

use it as it stands in subsequent calculations. 
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11.0 FINAL COMPUTATION 


11.1 Incompressible Option 


Once the velocity has been computed as described in Section 10.1, the only 
remaining calculation is the pressure coefficient 



Quantities output are described in Section 14.6. 


(56) 


11.2 Basic Compressible Calculation 


Once the calculations of Section 10.5 have been completed the following 
quantities are known at each on-body point, off-body point and control station 
point: 


Compressible Velocity Magnitude: 
Incompressible Velocity Vector: 
with components: 
and magnitude: 

Flow Passage Area: 

Average Incompressible Velocity: 

Average Density Ratio: 

These can now be used to compute other quantities 
Compressible Velocity Components: 

v y = v 1 y( v/v t ). 

Mach Number 


v i 

V V v iz 


App 

F i 

e = (p/p t ) 



(v/v.) 



(57) 


(58) 


(59) 
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Pressure Ratio 



( 60 ) 


The quantities output at all on-body, off-body, and control station points 
are described in Section 14.7. 


11.3 Additional Output Quantities Pertaining to a Cross Section as a Whole 

Additional overall quantities calculated are: 

Mass Flow Rate: 


w 


int 



(61) 


where ^ is compressible velocity with components computed in Section 11.2 
above, and the density ratio to be used in equation (6a) is 



Corrected Mass Flow Rate: 


(62) 


w 


_ *int^ 


cor 


(63) 


Specific Corrected Mass Flow Rate: 



ft 

A 


cor 

FP 


Average Mach number 


M is obtained by iteration using 

,3 


M = ^*- 


1 + 


(1/5)^ 


85.38 


with initial value M = 0. 


(64) 


(65) 
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11.4 An Example Comparing Calculated Results with Experimental Data 

To date only one realistic three-dimensional case has been calculated and 
the results compared with experimental data. The inlet is the scoop inlet shown 
in the sketch of Figure 11. Calculated and experimental surface pressures are 
compared along three curves at different "circumferential" locations around the 
inlet. Since the freestream is at zero angle of attack and yaw, if the inlet 
were axi symmetric, all three pressure distributions would be identical. Thus 
the variation between the three pressure distributions is a measure of the 
three-dimensionality of the flow. This variation is predicted rather well by 
the calculations, which also agree quite well with the data. Obviously, a good 
many more cases must be run to establish the accuracy of the present method, 
but the results shown in Figure 11 are certainly encouraging. 
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Figure 11. Comparison of calculated and experimental pressure distributions inside a scoop inlet. 


12.0 INPUT TO THE PROGRAM 


12.1 Types of Execution 

There are three modes of operation of the program: 

1. Complete Run (Mode 0) 

In this mode the program calculates the fundamental solutions, and 
continues to produce combined flow for prescribed values of the 
physical parameters. Input consists of geometric specification of 
the inlet including cross sections and a control station, and the 
inlet operating conditions. The fundamental solutions may be stored 
for future use as in Mode 1 below. It is anticipated that usually 
interrupted runs will be desired and that this mode will be run 
infrequently. 

2. Fundamental Solution Run (Mode 1) 

In this mode only the fundamental solutions are calculated and stored 
for future use on external units 4, 22, 23, 29, 32, 33 (see Appendix B). 
Cross sections, a control station, and operating conditions are not required 
This is the time-consuming portion of the calculation and is performed 
only once for each geometry (see Figure 5, Section 5.0). 

3. Combination Program Run (Mode 2) 

This is the most common type of run. Obviously either Mode 0 or Mode 1 
must have been run previously to produce the Fundamental Solution Data 
Set. This run starts with this data set and produces combined flow 
solutions. Input consists of cross sections, a control station, and 
the inlet operating conditions (see Figure 6, Section 5.0). 


12.2 The Various Input Cards 

The following are the input cards required by the program: 
Card #1: Option card 

Card #2: Title card 
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Card #3: 
Card #4: 
Card #5: 
Card #6: 
Card #7: 
Card #8: 
Card #9: 
Card #10 
Card #11 
Card #12 
Card #13 
Card #14 
Card #15 
Card #16 


Control data card 

Inlet strips control 

Vorticity coefficient card (optional) 

On-body coordinates card 

Off-body point card (optional) 

Viewing angle card 
Viewing angle input data 
Total cross-section point input 
Cross-section point data card 
Incompressible combination card 
Compressible combination Input card #1 
Compressible combination Input card #2 
Contour plot card (optional) 

Streamline tracing input card (optional) 


Obviously either Card 12 or Cards 13 and 14 are used. Cards from 10 on 
are unnecessary for Mode 1, and Cards 2 through 9 are unnecessary for Mode 2. 
However, Cards 10 and 11 and either 12 or 13 and 14 are necessary for every 
Mode 2 run. 


12.3 Specific Input Card Formats 


Card #1 — Option Card (Generally a nonzero integer calls for the option) 
Format (1315) 

Card 

Columns Variables 


1-5 

IFUND 

-It. II 

o o 

means do not calculate fundamental solutions, 
means do calculate fundamental solutions. The first 
is used only for Mode 2 execution 

6-10 

ICR0S 

o o 

II *1^ 

means no cross sections are input, 
means cross sections are input 

11-15 

ICMPRS 

-N 

where N is the number of compressible combinations 
to be run 

16-20 

IC0MB 

=M 

where M is the number of incompressible 
combinations to be run. 



Note: Either ICMPRS = 0 or IC0MB = 0 

For Mode 1, both are 0. 
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Card 

Columns Variables 

21-25 NCS = the number of the cross section to be used as the 

control station 

26-30 IPRINT =0 will print the combination program output 

=1 will suppress combination and cross-section output 
=4 will print cross-section geometry as well as 
combination output 

31-35 ISYMPR =0 means print solution only on input half of the inlet 

to means print solution also on the reflected half of 
the inlet (used with cases having yaw) 

36-40 INUNIT Unit that will read the input data (=5, if input is from 

cards) 

41-45 INSAVE =0 means fundamental solutions are not saved 

tO means save fundamental solutions and all other data 
necessary to restart with Mode 2 

46-50 NSL = number of streamlines to be used in the streamline 

tracing option 

51-55 IPERS =0 means use orthographic projection to plot streamlines 

7*0 means use perspective projection 

56-60 IS0PLT =0 means no cross-section plots 

7*0 means plot contours of constant pressure, etc. at 

each cross section 

61-65 IGE0M =0 means use input points unaltered 

7*0 means use geometry package to repanel 

Card #2 — TITLE card (Required if IFUNDj*0 on card #1) 

Format (18A4) 

Card 

Col umn Variables 

1-68 TITLE(l) An alphanumeric array which describes the case data to 

to be run. The array should contain no more than 68 

TITLE(17) characters. 

69-72 TITLE(18) = '0TW0' means two points per card format is used for 

inputting body coordinates 
= blanks means one point per card format is used 
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Card #3 — Control Data Card (Required if IFUND^O on card #1 ) 
Format (A4, 3X, 213, 12X, 513, 6X, 4L1, 2F10.0) 

Card 


Columns 

Variables 

1-4 

CASE 

Case ID number 

8-10 

N0FF 

=0 means no off-body point is input 
i*0 means there is off-body point input 

11-13 

LIST 

=0 program continues after body panels are formed 
j*0 program terminates after body panels are formed 
(i.e. program ends at the completion of INPUT 
subroutine) 

26-28 

NSYM 

=0 symmetry is not used, entire body is input 
=1 one plane of symmetry is used. Only half the 
body, for which y ^0, is input. 

29-31 

IB 

Input vorticity coefficient option 

=0 no B-values input. Program sets all B's to 1.0 

=1 uses input B's. 

32-34 

N0RM 

=1 lengths will be normalized by an input reference 
length 

=2 velocities will be normalized by input reference 
velocity 

=3 both lengths and velocities will be normalized 

35-37 

METIN 

=1 means the input is in metric units 
=0 means the input is in English units 

38-40 

MET0UT 

=1 means metric output is desired 
=0 means the output is in English units 

47 

TYPEA 

=T if there is no centerbody 
=F if there is a centerbody 

48 

TYPEB 

=T if the front of the centerbody is aligned with 
the inlet highlight 
=F otherwise 

49 

TYPEC 

=T if the centerbody is completely iiu ;de the inlet 
downstream of the highlight 
=F otherwise 

50 

TYPED 

=T if the centerbody protrudes outside the inlet 
upstream of the highlight 
=F otherwise 


(Note: 

Exactly one of the previous four flags must be set to T, 
the other three to F.) 
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Card 

Columns Variables 


51-60 

SF 

scale factor 

61-70 

RFLENG 

reference length for normalizing the length data. 
Input only if N0RM=1 or N0RM=3 

Card #4 — Inlet Strips Control Card (Required if IFUND^O on card #1) 
Format (314) 

Card 

Columns Variables 

1-4 

NS0RCE 

number of on-body source panels in each strip of the 
inlet (lifting body) section 

5-8 

NWAKE 

number of wake panels in each strip of the inlet. 
Usually one 

9-12 

NSTRIP 

number of strips of panels on the inlet 


Card #5 — Vorticity Coefficient Input Card (This card is input only if IB=1 is 
input on card #3 and IFUND^O on card #1.) 

Format (6 FI 2. 6) 

Card 

Col umns Variables 

1-12 B(I) values of vorticity coefficient —6 values per card. 

13-24 Repeat the same format on next card if more than 6B's 

25-36 are input. 

37-48 (Note: The total number of B's input must equal the 

49-60 total number of strips that define the inlet.) 

61-72 

Card #6 - On-Body Coordinates Input Card (Required if IFUND^O on card #1.) 

Format (3E10.0, 211) — For one point per card 
Format (2(3E10.0, 211)) - For two points per card 

The maximum number of on-body panels on inlet and centerbody together 
is 1000. The maximum number of input points may slightly exceed this. 


Card 


Columns 

Variables 

1-10 

X(I) 

11-20 

Y(I) 

21-30 

Z(I) 


The x,y,z coordinates of an input point on an N-line, 
either on-body or wake (inches in English units) 
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Card 

Columns 


Variables 


31 ISTAT Status flag 

=2 new section (1st point of inlet and of centerbody) 
=1 new N-line 
=0 same N-line 
=3 end of all input data 

(Note: The status flags are somewhat different if the geometry 

package is used to repanel the body, Section 13.0.) 

32 LABEL =1 lifting body - used for the inlet 

=0 nonlifting body — used for the centerbody 

(Note: For two points per card repeat the above sequence from 

Card Columns 33 - Card Columns 64.) 

Card #7 — Off-Body Point Cards (Required if N0FF? { O on card #2 and IFUND^O on 

card #1) 

Format (2(3E10.0, II)) 

Card 

Col umns Variables 


1-10 

X0FF ( I ) 


of an input off-body point 

11-20 

Y0FF(I) 

x,y,z coordinates 

21-30 

Z0FF(I ) 



31 

ISTAT 

Status flag 
=3 indicates end 

of all off -body point input 


(Note: 

Repeat the above sequence from Card Columns 32 - Card 


Columns 63. ) 

Card #8 — Total Number of Viewing Angles Input Card (Required if IFUND^O) 
Format (15) 

Card 

Columns Variables 

1-5 NVIEW number of viewing angles input 


Card #9 — Viewing Angle Input Data Card (Required if IFUND^O) 

Format (3F6.0, 12, 15A4) One viewing angle is input per card (required even if 
graphics is not called for) 


60 



Card 


Columns 

Variables 


1-6 

7-12 

13-18 

PSI 

THETA 

PHI 

angles the viewing direction makes with the x,y,z 
directions, respectively (degrees) (see Section 15.1) 

19-20 

LAST 

set equal to 1 if this is the last viewing angle card, 
otherwise, set it to 0 

21-80 

PLABEL 

a 60-character array that describes the viewing angle 

Card #10 ■ 

- Total Number 

of Cross-Sections Input Card (Required if ICR0S#O) 

Format (15) 


Card 

Columns 

Variables 


1-5 

NPLANE 

number of cross-section points input 

Card #11 • 

- Cross-Section Input Data Card (Required if ICR0S#O) 

Format (6E10.0, 15) 


Card 

Columns 

Variables 


1-10 

XCSP 


11-20 

YCSP 

x,y,z coordinates of the input cross-section point 

21-30 

ZCSP 

(inches in English units) 

31-40 

UI 


41-50 

UJ 

components of the unit vector of the input cross-section 

51-60 

UK 

point. 

61-65 

ICSR 

number of panels In each of the strips in the cross- 
section plane (a maximum of ten) 

Repeat this format for 

another cross-section point until the total number of 

cards #11 

equals NPLANE 

input in card #10. 

Card #12 

- Combination 

(Incompressible) Program Input Card (This card is required 


if IC0MB^O on card #1.) 

Format (5F10.6, 215) 


Card 

Columns 

Variables 


1-10 

ALFANG 

angle of attack a in degrees 

11-20 

YAWANG 

yaw angle B in degrees 
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Card 

Columns 


Variables 


21-30 VINF freestream velocity V 

31-40 VREF reference velocity 

41-50 VAVG average velocity V 

V 

51-55 ISKIPC =0 combination program will be called 

=1 bypass the combination program even though IC0MB? f O, 
and go directly into the streamline tracing program 

56-60 IDEL =0 use the program default delta values for cross- 

section plot 

=1 input delta values to be used for cross-section plot 
Card #13 — Flow Passage Program Input Cards (Required only if ICMPRS^O on card #1.) 

(a) Format (15) 

Card 

Columns Variables 

1-5 NCPT number of centerline points input 

(b) Format (6E10.0) 


(feet per second 
in English units) 


Card 

Columns Variables 


1-10 XCPT 
.11-20 YCPT 
21-30 ZCPT 


x,y,z coordinates of the centerline point 
(inches in English units) 


31-40 UI 
41-50 UJ 
51-60 UK 


components of the unit normal vector of the 
centerline point 


Repeat card (b) as many times as specified by NCPT in card (a). 


The following set of cards is also needed if the inlet is TYPE D, i.e. the 
centerbody is upstream of the inlet lip. 

(c) Format (15) 

Card 

Columns Variables 

•1-5 NE6CPT number of centerline points upstream (outside) of 

the inlet lip 
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(d) Format (6E10.0) 


Card 



Col umns 

Variables 


1-10 

XNEG 


11-20 

YNEG 

x,y,z coordinates of an upstream centerline point 

21-30 

ZNEG 


31-40 

VI 


41-50 

. VJ 

components of the unit vector of an upstream centerline 

51-60 

VK 

point 


Repeat card (d) as many times as specified by NEGCPT in card (c). 

Card #1 4a — Combination (Compressible) Program Input Card (This card is required 
if ICMPRS#0 on card #1 . the format and input variables are the same 
as in card #12.) 

Format (5F10.0, 215) 

Card 

Col umns Variables 

1-10 ALFANG angle of attack a in degrees 


11-20 YAWANG 

21-30 VINF 

31-40 VREF 

41-50 VAVG 

51-55 ISKIPC 

56-60 IDEL 


yaw angle e in degrees 
freestream velocity V ot 
reference velocity V ^ 
average velocity V c 

same meaning as in card #12 


(feet per second 
in English units) 


Card #1 4b - Second Input Card for the .Compressible Combination Program 


Format (7F10.0) 


Card 

Columns 

Variables 


1-10 

RMINF 

Mach number M 

CO 

11-20 

RMAVG 

average Mach number M 

21-30 

PT0TAL 

total pressure P^ 

31-40 

PSTATC 

freestream static pressure 



see Section 10.0 
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Card 

Columns 


Variables 


41-50 

51-60 

61-70 


TEMPT 

TEMPS 

FRATE 


total temperature T^ 
freestream static temperature T g j 
inlet mass flow rate w 


(degrees Rankine) 
in English I 
uni ts ) > 


lb mass/sec 


} See 
i Secti 
) 10.0 


on 


Card #15 - Increment Values for Contour Plots (1 card for each input cross-section 
point.) (Required if and only if IDEL=1 ) 

Format (4F10.0) 


Card 

Columns 

Variables 


1-10 

DELCP 

isobar increment (default value = 0.1) 

11-20 

DELTHX 

flow angle in the r.-di recti on contour plot increment 
(default value = 1°) 

21-30 

DELTHY 

flow angle in the n-di recti on contour plot increment 
(default value = 1°) 

31-40 

DELMCH 

Mach number increment (omit this for incompressible 
combination) (default value = 0.1) 

Card #16 - 

Streamline 

Tracing Input Data Card(s) (One card per streamline) 

Format (415 

, F10.0, 215, 311, 2X, 311) 

Card 
Col umns 

Variables 


1-5 

ICFLAG 

centerbody or inlet flag: 1 -*• centerbody; 2 inlet* 

6-10 

IPIESL 

strip number on the centerbody (or inlet) 

11-15 

IELN0 

the panel number of this strip** 

16-20 

IDIREC 

the direction of the streamline march: +1 -*■ in local 

flow direction; -1 -*■ opposite to local flow direction 

21-30 

XCUT0F 

maximum x-value for streamline calculation (cutoff 
condition) 

31-35 

NPMAX 

maximum number of streamline points to be allowed for 
this streamline (cutoff condition) 


♦centerbody (if any) must be input before the inlet if streamlines are to be 
calculated. 

**streamline starts at the control point of this panel. 
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Card 


Columns 

Variables 


36-40 

IPRINT 

print flag for streamline tracing calculations: 

0= minimum print. (This is the normal option for 
production case) 

2= maximum print (for batch or TSO debug) 

3= medium print (mostly for TSO debug) 

41 

IMVSS 

Flag to generate a plot of Mach No. vs normalized 
distance along this streamline. 

42 

IPVSS 

Flag to generate a plot of Ps/Pt vs normalized 
distance along the streamline. 

43 

IAVSS 

Flag to generate a plot of flow inclination angles 
(rel. to x-axis) vs normalized distance along this 
streamline. 

44-45 

(blank) 


46 

IMVSX 

Flag to generate a plot of Mach No. vs x for this 
streamline. 

47 

IPVSX 

Flag to generate a plot of P s /Pt vs x for this 
streamline 

48 

IAVSX 

Flag to generate a plot of flow inclination angles 
(rel. to x-axis) vs x for this streamline. 


Card no. 16 should be repeated NSL times. 
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13.0 INPUT TO THE GEOMETRY PACKAGE 


Note: If it is desired to use the geometry package (references 6 and 7) 

to repanel the body, then IGE0M^O is specified on card #1. The 
following instructions pertain to the input cards for the geometry 
package. The geometry package has been greatly simplified for use 
in the inlet program, but it is still more general than is needed. 
For example, many of the spacing options are superfluous. The 
inlet and centerbody are referred to as body 1 and body 2. Either 
may be input first*. They should be inserted behind card #5 and 
card group #6 is not needed if this option is chosen. 


Required 

Card or Card 


Type 

Optional 

Col umns 

Format 

Variable 

Description 

1 

Requi red 

1-4 

14 

NSEC 

Total number of sections 

2 

Required 

1-4 

14 

NB(I ) 

Specified number of on-body 
points per N-line on body I 



5-8 

14 

Ntf(I) 

Specified number of wake points 
per N-line on body I 



9-12 

14 

NS (I ) 

Specified number of N-lines on 
body I 



13-16 

14 

NALGB (I) 

Spacing algorithm for on-body 
points on N-lines 


NALGB=0, Input distribution, 
unaltered 

NALGB=1 , Input distribution, 
augmented 

NALGB=2, Constant increments 
in arc length 

NALGB=3, Cosine distribution 

NALGB=4, Curvature-dependent 
distribution 

NALGB=5, User-specified dis- 
tribution 

17-20 14 NALGW(I) Spacing algorithm for wake points 

on N-lines. (Values have same 
significance as values of NALGB, 
except that NALGW=3 and NALGW=4 
should not be used.) Normally 
only one wake point is used. 

♦However, if streamlines are to be calculated, the nonlifting body must be input 
fi rst. 
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Required 

Card or Card 

Type Optional Columns Format Variable 


2 Required 21-24 
(cont) 


14 NALGS(I) 


Description 


Spacing algorithm for N-lines 
NALGS=0, Input distribution, 
unaltered 

NALGS=1 , Input distribution, 
augmented 

NALGS=2, Constant increments 
NALGS=3, User-specified dis- 
tribution 


29-32 

14 

NSEG(I) 

Number of segments into which 
M-lines are broken (Default=l) 
This is to allow for corners 

33-36 

14 

NTR(I) 

Component transformation flag. 
(The value indicates the number 
of transformations to be per- 
formed.) 

37-40 

14 

IP1(I) 

Punch flag for transformed input 
coordinates 
IP1=0, No punch 
IPl^O, Punch 

41-44 

14 

IP2 (I ) 

Punch flag for coordinates after 
repaneling bodies 
IP2=0, No punch 
IP2^0, Punch 

4 Optional 1-10 

11-20 
21-30 
31-40 
41-50 
51-60 
61-70 
71-80 

8F10.0 

SSB ( 1 ,1 ) 
SSB ( I ,2) 

* 

• 

Specified arc lengths of on-body 
points on N-lines on body I 
0.0 £ SSB ^ 1.0 
8 points per card, NB(I) total 

This card is required if 
NALGB ( I ) =5 
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Required 

Card 




Card 

or 



Description 

Type 

Optional 

Columns 

Format 

Variable 

5 

Optional 

1-10 

8F10.0 

SSKCI.U 

Specified arc lengths of wake 


11-20 


SSW(.I,2) 

points on N-lines on body I 



21-30 


• 

0.0 < SSW < 1.0 



31-40 


• 

8 points per card, NW ( I ) total 



41-50 

51-60 

61-70 

71-80 


• 

This card is required if 
NALGW(I)=5 

6 

Optional 

1-10 

8F10.0 

SSS(I,1) 

Specified distribution of N-lines 


11-20 


SSS(1 ,2) 

on body I 



21-30 


• 

0.0 < SSS < 1.0 



31-40 


• 

8 points per card, NS (I ) total 



41-50 

51-60 

61-70 

71-80 


• 

This card is required if 
NALGS(I )=3 

7 

Optional 

1-4 

514 

IEND(I,1) Point numbers of the ends of the 


5-8 


IEND(I,2) segments on M-lines on body I 



9-12 


• 

NSEG(I) total (max. =5 



13-16 


• 




17-20 


• 

This card is required if 
NSEG ( I ) >1 

8 

Optional 

1-4 

14 

ITR2 (I ,J) 

Type of Jth transformation to be 
applied to body I 






ITR2=1 , Scaling 
ITR2=2, Translation 






ITR3=3, Rotation 



11-20 

F10.0 

TR(I,J,1) 

First transformation parameter: 


x multiplication factor if ITR2=1 
x translation if ITR2=2 

angle of rotation if ITR2=3 

(positive clockwise, looking down 
axis cf rotation away from origin) 
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Card 

T ype. 

8 

(cont) 


Required 

or Card 

Optional Columns' ~Fontiat Variable- Description 

Optional 21-30 F10.0 TR(I,J,2) Second transformation parameter: 

y multiplication factor if ITR2=1 
y translation if ITR2=2 

x direction cosine of axis 
of rotation if ITR2=3 

31-40 FI 0.0 TR(I,J,3) Third transformation parameter: 

z multiplication factor if ITR2=1 
z translation if ITR2=2 

y direction cosine of axis 
of rotation if ITR2=3 

41-50 FI 0.0 TR(I,J,4) Fourth transformation parameter: 

Dummy if ITR2=1 or ITR2=2 
z direction cosine of axis of 
rotation if ITR2=3 

Note : The axis of rotation is 

assumed to pass through the 
origin. Up to 10 transformations 
may be performed, in any order. 
Input this card type NTR(I) times. 


Required 1-10 

F10.0 

X 

11-20 

FI 0.0 

Y 

21-30 

F10.0 

Z 

31 

11 

STATUS 

32 

11 

LABEL 


x,y,z coordinates of a point (all 
points input N-line by N-line, 
body by body). 

STATUS=0 same N-line 
STATUS=1 new N-line 
STATUS=2 new body 
STATUS=3 last point input (excep- 
tion is #5 below) 
STATUS=4 first point on wake 
STATUS=5 first point on wake and 
also last point input 

LABEL=0 nonlifting body 
(center body) 

LABEL- 1 lifting body (inlet) 
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33-42 

F10.0 

X 

43-52 

F10.0 

X 

53-62 

FI 0.0 

1 

63 

11 

STATUS 

64 

11 

LABEL 


x,y,z coordinates of the next 
point input 


Same meaning as above 


Repeat cards to inout more 
points. If TITLE(18HTW0 on 
card no. 1, input 2 points per 
card. If TITLE(18) is not equal 
to "0TW0," it implies 1 point 
per card. 



14.0 TABULAR OUTPUT 


14.1 Basic Body Geometry 

N Strip number 

M Panel number 

X, Y, Z Coordinates of the four input points used to form the panel 

NX, NY, NZ x,y,z components of the panel normal vector 
X0, Y0, Z0 Coordinates of the control point 

D For lifting (inlet) panels, it is an average panel length 

For nonlifting (centerbody) panels, it is the projection 
distance required to make the input points coplaner 

T Maximum diagonal of the panel 

A Area of the panel 

14.2 Fundamental Solution Output 

I Panel number. Here all panels are listed in a single array 

containing all strips 

STATVX x,y,z components of the velocity at a control point for the 

STATVY static solution 

STATVZ 

UVX1 x,y,z components of the velocity at a control point for the 

UVY1 (1,0,0) flow, in which the uniform freestream is parallel to 

UVZl the x-axis 

UVX2 x,y,z components of the velocity at a control point for the 

UVY2 (0,0,1) flow, in which the uniform freestream is parallel to 

UVZ2 the z-axis 

UVX3 x,y,z components of the velocity at a control point for the 

UVY3 (0,1,0) flow, in which the uniform freestream is parallel to 

UVZ3 the y-axis 
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14.3 Cross-Section Control Data 


C.S.N0. 

XCSP 

YCSP 

ZCSP 

ur 

UJ 

UK 

RADIAL 

PANELS 


J 

I 

X, Y, Z 

XC 

YC 

ZC 

XI 

ETA 

AREA 


Cross-section number 

Input coordinates of the cross-section point 

Input components of the cross-section normal vector 

Number of panels to be formed on each radial strip in the 
cross-section plane 

14.4 Cross-Section Panel Geometry 

Strip number 

Panel number 

Coordinates of the corner points of a panel 

x,y,z coordinates of the panel control point where flow 
quantities are to be calculated 

Coordinates of the control point in the cross-section 
coordinate system 

Panel area 
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14.5 Cross-Section Flux 


FL0W Indicates what type of flow 

FL0W=1 static solution 

FL0W=2 (1,0,0) flow, uniform freestream parallel to the x-axis 

FL0W=3 (0,0,1) flow, uniform freestream parallel to the z-axis 

FL0W=4 (0,1,0) flow, uniform freestream parallel to the y-axis 

FLUXX Components of the vector volume flux at the cross-section, 

FLUXY equations (28) and (29) plus subsequent discussion 

FLUXZ 

T0TAL FLUX Normal volume flux through the cross-section plane, equation (29) 

14.6 Incompressible Combination Output 
N Strip number 

M Panel number 

XC 

YC Coordinates of control point 

ZC 

NX 

NY Components of normal vector 

NZ 

C0MBVX 

C0MBVY Components of velocity 

C0MBVZ 

VN Normal velocity 

V Total velocity magnitude 

CP Pressure coefficient 

DCX 

DCY Components of directional cosines of velocity 

DCZ 

The above explanations applied to all elements (i.e. on-body, off-body & 
cross-sections) . 
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14.7 Compressible Combination Output - On Body 
N Strip number* 

M Panel number 

XC, YC, ZC Control point coordinates 

VX, VY, VZ Velocity components 

VMAG Velocity magnitude “vVy +V** +V^ 

x y z 

MACH Local Mach number 

P/PT Static-to-total pressure ratio 

VBAR Average incompressible velocity 

EPSL0N e, Average density ratio, equations(53) and (54) 

AREA (FP) Flow passage area 

*If letter R is suffixed to the strip number, it denotes that the strip is 
on the reflected half of the body. 

14.8 Error Messages 

(1) Message: **END OF INPUT DATA. NORMAL PROGRAM TERMINATION. 

Source: MAIN program 

Cause: End of all data. 

Action: None 

(2) Message: THE MATRIX OF XXXXX ELEMENTS EXCEED XXXXXX. 

Source: MAIN program 

Cause: The input matrix is too large for the program's allowable 

computer storage. 

Action: Change the working area size in MAIN program. 

(3) Message: TOTAL NUMBER OF CENTER-BODY INTERSECTION POINTS DOES NOT 

EQUAL TO THE NUMBER OF NONLIFTING NLINES INPUT. 

Source: Subroutine CSGE0M 

Cause: Each cross-section plane (defined by the input cross-section 

point and its normal vector) should give one intersection 
point with each line that defines the centerbody. If the total 
number of intersection points does not equal the number of 
input nonlifting N-lines, the above message will be printed and 
the program continues. 



Action: Check the input cross-section point and its normal vector to 

see if the cross-section plane does cut through each nonlift- 
ing N-line. 

(4) Message: (a) TOTAL NUMBER OF INTERIOR POINTS DOES NOT EQUAL TO THE 

LIFTING NLINES INPUT. 

(b) TOTAL NUMBER OF EXTERIOR POINTS DOES NOT EQUAL TO THE 
LIFTING NLINES INPUT. 

Source: Subroutine CSGE0M 

Cause: Intersection of the cross-section plane and the lifting N-line 

should generate two intersecting points as illustrated below: 


LIFTING 

N-LINE 


If the number of the interior point or the exterior point 
does not equal the number of lifting H-lines, the above 
message will be printed but program continues to execute. 

Action: Check the input cross-section point and normal vector to see 

if the cross-section plane does cut through the lifting N-line. 

(5) Message: CURRENTLY EXECUTING SECTION TYPE = , XXX 

SECTION TYPE FROM UNIT NO. 4 =, XXX 
PROGRAM STOPS 

Source: Subroutine FLPASG 

Cause: External unit No. 4 stores the geometrical quantities of the 

input panels for both the inlet and centerbody. If the body 
type retrieved from unit 4 does not match the body type that 
the program is currently using, the above message will be 
printed and execution stops. 

Action: Check the contents of the external unit No. 4 and see if it 

matches the case that is used. 
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(6) Message: NON-LIFTING INTERSECTION POINTS NOT EQUAL TO THE NUMBER OF 

NLINES INPUT 

Source: Subroutine FLPASG 

Cause: 1 

> See Message No. (3) 

Action: ) 

(7) Message: INTERIOR POINTS NOT EQUAL TO THE NUMBER OF NLINES INPUT. 

Source: Subroutine FLPASG 

Cause: ) 

> See Message No. (4) 

Action: ) 

(8) Message: EXTERIOR POINTS NOT EQUAL TO THE NUMBER OF NLINES INPUT. 

Source: Subroutine FLPASG 

Cause: } 

| See Message No. (4) 

Action: ) 

(9) Message: ALLOWABLE NUMBER OF POINTS EXCEEDED IN SUBROUTINE PANELS, 

EXECUTION TERMINATING. 

Source: Subroutine GEINPT 

Cause: Input total number of panels (for both inlet and centerbody) 

exceeds 1100. 

Action: Change the input geometry and rerun. 

(10) Message: MISMATCH OF ELEMENTS IN A LIFTING STRIP IS DETECTED. 

ELEMENTS FORMED = XXX, ELEMENTS INPUT = XXX. 

COMPUTATION TERMINATED. 

Source: Subroutine INPUT 

Cause: Input number of strip and elements for the inlet is not con- 

sistent with the input lifting information card. 

Action: Check the body coordinates input cards and see if the flags 

for new line or new section are punched correctly. 

(11) Message: NPMAX=XXX HAS BEEN REACHED. COMPUTATIONS TERMINATING FOR THIS 

STREAMLINE. 

Source: Subroutine TRACE 
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Cause: Generated streamline tracing exceeds the allowable dimension. 

Program goes to the next streamline. 


Action: None 


(12) Message: INPUT XX PLANE OF SYMMETRY OPTION IS NOT IN THE PROGRAM. 

IT IS REPLACED BY 1 PL. SYM. 


Source: 
Cause: 
Action: 
(13) Message: 


Subroutine VFMLFT and VFMNLF 

Input number of symmetry plane is more than 1. 

User's action none. Program changes the input symmetry to one. 

ERROR IN VFORM. THE ELEMENTS FORMED DO NOT CORRESPOND TO THE 
NO. OF BODY ELEMENTS. 


Source: Subroutine VFMLFT and VFMNLF 


Cause: 
Acti on : 

(14) Message: 
Source: 
Cause: 

Acti on : 

(15) Message: 
Source: 
Cause: 
Action: 

(16) Message: 

Source: 

Cause: 


On-body panel counter error 

Check and see if the variables appeared in the labelled 
commons "INFORM" and "INLIFT" have been changed. 

LABEL ERROR IN LIFTING VFORM 

Subroutine VFMLFT 

Data on the external unit No. 4 did not correspond to the type 
of body at the time of execution. 

Check the data on unit No. 4. 

LABEL ERROR IN NONLIFTING VFORM. 

Subroutine VFMNLF 

Same as in (14) 

Same as in (14) 

TOTAL NUMBER OF CONTROL POINTS SHOULD BE XXXXX. 

THE NUMBER OF POINTS STORED ON UNIT NSETV IS XXXXX. 

EXECUTION STOPS. 

Subroutine C0MPR0, C0MPRS 

The total number of control points (on-body plus off-body) in 
the variable N0N is not equal to the variable NP0INT stored in 
the external unit NSETV. (NSETV stores all of the fundamental 
solutions. ) 


77 



Action: 

(17) Message: 

Source: 

Cause: 

Acti on : 

(18) Message: 
Source: 
Cause: 


Acti on : 
(19) Message: 


Source: 

Cause: 


Action: 

(20) Message: 

Source: 

Cause: 


Check and see if the unit NSETV did store the correct data. 

INPUT CONTROL STATION NO. = XXXXX. 

READ FROM UNIT NFLUX STATION NO. = XXXXX. 

PROGRAM STOPS FOR CHECKING. 

Subroutine COMPRO, COMPRES 

Input control station number NCS is not the same as the control 
station number stored in unit NFLUX. 

Check contents of unit NFLUX. 

ON-BODY ELEMENT MISSED. JOB STOPS. 

Subroutine COMPRO, COMPRS 

Total number of on-body control panels is stored in variable 
KONTRL. The counter KK is used for counting total number of 
on-body panels. When these two variables are not equal, the 
message will be printed. 

Check labelled common CONFLG to see if the variable KONTRL 
has been changed. 

CROSS-SECTION PLANE NUMBER CURRENTLY USING XXXXX PLANE NO. 

FROM UNIT NV IS XXXX. 

PLANE NO. FROM UNIT NOFSAV IS XXXX. 

STOP 77. 

Subroutine COMPRO, COMPRS 

The three numbers explained in the message should be the same, 
if they are not, it indicates either the units have different 
contents or the counter for the cross-section plane number IP 
is incorrect. Stop code 77, 78 applied to COMPRO and COMPRS, 
respectively. 

Check external units NV and NOFSAV. 

STRIP NO. NOW USING IS XXXX. 

STRIP NO. FROM UNIT NOFSAV IS XXXX. 

PROGRAM HALTS. 

Subroutine COMPRO, COMPRS 

Unit NOFSAV is used in retrieving the cross-section geometry 
data. Strip counter JI is read from NOFSAV and do- loop (2800 
in COMPRO, 3800 in COMPRS) counter J is the current strip 
counter. If these two variables (Jl and J) are not equal, it 
indicates data read in from NOFSAV is incorrect. 
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Action: Check contents of the unit NOFSAV. 

(21) Message: POINT NUMBER ERROR. ON UNIT NV N = XXXXX. 

BUT N IS INCREMENTED = XXXXX. 

STOP 59. 

Source: Subroutine COMPRO, COMPRS 

Cause: In do-loop (2700 in COMPRO, 3700 in COMPRS) counter JS is 

incremented for each control panel. Unit NV also has a 
counter N which should be the same as JS. If these two do 
not agree, message is printed. 

Action: Check contents of unit NV. 

(22) Message: ********** EPSLON AT 100TH ITERATION IS USED **********. 

Source: Subroutine COMPRS 

Cause: Iteration for the computation of e equal to 100. The 100th 

iterated e value is used. 

Action: None 

(23) Message: ********** AVERAGE MACH NO. AT 100TH ITERATION IS USED ********** 

Source: Subroutine COMPRS 

Cause: Iteration for the computation^ of the average Mach number M equal 

to 100. The 100th iterated M is used. 

Acti on : None 

(24) Message: AFTER 100 ITERATIONS, COMPUTED VC IS +XXXXXX.XXXXXXXXX DELTA 

VC IS ±XXXXXX.XXXXXXXX. 

PROGRAM STOPS. 

Source: Subroutine COMPRS 

Cause: Computed SV C > 10“^, program stops. 

Action: Reconsider the input data to the combination program (compres- 

sible) and make proper change. 

(25) Message: NUMBER OF (PRESSURE MACH NO. FLOW ANGLE) ISOLINES EXCEEDS 

DIMENSION ALLOWED. PROGRAM STOPS. 

Source: Subroutine ISOPLN. 

Cause: Computed total number of isolines exceeds the dimension allowed 

in the program. 

Action: Change the input delta values for the isoline plot and reinput. 
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15.0 EXAMPLES OF GRAPHICAL OUTPUT 


15.1 Surface Streamlines 

The user identifies a particular body panel by inputting the number of 
the strip on which it lies and the number of the panel on that strip. The panel 
may be on either the inlet or the centerbody. The program then calculates the 
surface streamline that goes through the control point of that panel in either 
the upstream or downstream direction. (If both are desired, two inputs are 
required, since the program regards them as two separate streamlines.) A number 
of panels may be specified and streamlines calculated through each of their 
centroids. 

To facilitate interpretation of the results, two types of graphical output 
have been provided. One is rather straightforward. Local values of pressure, 
Mach number, and flow angle are displayed as functions of arc length (and/or 
Cartesian x) along the streamline on conventional two-variable Cartesian plots. 

The second type of graphical output is more complicated. The program plots 
the actual streamlines themselves against the panel representation of the inlet 
configuration. The angles from which the resulting configuration is viewed may 
be specified by the user. Their choice is quite important. Viewed from an 
inappropriate direction, the streamline pattern can be quite deceptive. A set 
of default viewing angles have been built into the program. In complicated cases 
more than one run may be required to determine the best angles. One complica- 
tion is the fact that only the panels facing the observer are plotted, while the 
streamline is plotted regardless. Thus a view of an interior streamline viewed 
from the exterior is quite misleading. 

The plotted view of the inlet with its streamlines may be obtained either 
by orthographic (parallel) projection or by perspective projection for which a 
viewing distance has been built into the program. Examples of the output are 
presented in Figure 12, which shows streamlines on the interior surface of an 
inlet at angle of attack for both types of projection. The viewing angle is 
10 degrees with respect to the inlet centerline, which often has been found to 
be an effective angle from which to view. While the viewing angle is the same 
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(a) Orthographic projection. 


PANEL EDGES 
STREAMLINES 


(b) Perspective projection. 



Figure 12. Calculated streamlines on the interior surface of an inlet. 


for both projections, it appears to be quite different. In this case the 
perspective view appears more informative, but examples can be found for which 
orthographic is to be preferred. 

The viewing angles input by the user are: a yaw angle PSI, a pitch angle 

THETA, and a roll angle PHI. They are used both for examining the geometry 
before performing the flow calculation and for viewing final results as dis- 
cussed in this section. Since the rotation axes pertinent to the viewing angles 
remain fixed in the body, it is simpler to imagine that the body and the 
coordinate system remain fixed while the observer is rotated by the prescribed 
yaw, pitch, and roll angles about the coordinate axes. In each case rotation 
about a particular coordinate axes implies that the value of that particular 
coordinate of the observer remains constant, as does his distance from the axes, 
and the positive sense of rotation is clockwise. Initially, the observer is 
assumed to be situated on the positive x-axis and to be looking towards the 
origin, which he continues to do at all times. He is first rotated about the 
z-axis by the yaw angle and thus remains in the xy-plane. The observer is then 
rotated about the y-axis by the pitch angle and finally about the x-axis by the 
roll angle. This last rotation is rarely used, because the observer can be 
brought to any desired inclination using only the yaw and pitch angles. 
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15.2 Equi -Value Contours in a Cross Section 


As described in Sections 9.0 and 14.0, plane cross sections are defined at 
various locations inside the inlet by means of a point and a normal vector. 

Control points are automatically distributed "spiderweb style" over the cross 
section from the input point to the inner surface of the inlet if there is no 
centerbody at that location or from the surface of the centerbody to the inner 
surface of the inlet if there is a centerbody. Coordinates are defined in the 
cross-section plane having the input point as origin. Tabular output of 
various flow quantities at the control points is provided. 

The graphical output option for a cross section first plots the geometric 
shape of the cross section using its own coordinate system, i.e., the intersection 
of the cross-section plane with the centerbody (if any) and the interior surface 
of the inlet. It then constructs four separate plots against this geometry 
consisting of curves of equal value of: nondimensional pressure (pressure coef- 

ficient incompressibly and static- to- total pressure ratio compressibly), Mach 
number, and flow angles with respect to the cross-section normal vector in the 
planes defined by this vector and the cross-section coordinate axes. Increments 
of value to be used for these curves may be input by the user, but the same 

values must be used for all cross sections. Default values for the plotting 

increments are 0.1 for nondimensional pressure and Mach number and one degree for 
the two angles. As an aid to interpretation, the equi-value curve having the 
smallest value of the quantity in question and the curve having a value halfway 
between the largest and smallest values of the quantity are plotted as heavy and 
dotted curves, respectively. The values of the quantity corresponding to these 
two curves and the value increment are printed on the plot. 

To illustrate this capability, the example shown in Figure 13 was selected. 

It is a round inlet but was panelled as a three-dimensional body. The locations 

of the two cross sections where flow quantities are computed are shown. The 
forward cross section was designated the control station where the average 
velocity was specified as three times freestream. The angle of attack was 40°. 
Figure 14 shows isobar plots at the two cross sections. Figure 15 shows plots 
of equal values of flow angle. Here the cross-section normal vector is parallel 
to the x-axis. Thus the flow angle in the xz-plane has the meaning of a pitch 
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17.0 PRINCIPAL NOTATION 


area of a panel 

area of the intersection of a plane with the centerbody 

area of the intersection of the initial centerline plane with 
the interior of the inlet 

flow passage area 

geometric area associated with exterior points of the inlet's 
surface 

normal velocity at ith control point due to unit source density 
on jth panel 

speed of sound, a* is critical speed, equation (43) 

combination constant for static solution, equation (5) 

strength of bound vorticity on the kth strip of the inlet 

subscript used with V,y,p,q, and P to denote average condi- 
tions over a cross section 

vector volume flux through a cross section 

subscript denoting quantities associated with the jth panel 

number of lifting strips 

Mach number 

total number of panels on both inlet and centerbody 

components of unit normal vector 

pressure 

control station mass flow used in combination program, equations 
(5), (30), and (51) 

number of control points in radial direction across a cross 
section 

dynamic pressure 

as a subscript denotes freestream static conditions; as a super- 
script denotes quantities associated with solution to flow about 
inlet in static operation 



number of N-lines on the inlet; also denotes temperature 
a subscript denoting reservoir conditions 
fluid velocity 
freestream velocity 

velocity vector at ith control point. Subscript 0 denotes 
onset flow velocity. Subscript k denotes velocities due to 
vorticity solution on kth strip. Superscript f denotes 
velocity of fth fundamental flow solution. Specifically f 
can be either: s denoting static solution, (100) denoting 

solution at zero angle of attack and yaw, (010) denoting pure 
yaw, or (001) denoting 90° angle of attack with no yaw. 

Cartesian coordinates 

coordinates of input cross-section point 

coordinates of intersection of kth centerbody N-line with cross- 
section plane 

coordinates of intersection of kth inlet N-line with cross- 
section plane 

angle of attack 

angle of yaw 

density ratio, equation (53) 

Cartesian coordinates in the plane of a cross section. The 
n-axis is parallel to the y = 0 plane, which is the symmetry 
plane in symmetric cases. 

fluid density 

source density 



18.0 REFERENCES 


1. Stockman, N.O. : Potential and Viscous Flows in VTOL, STOL or CTOL 

Propulsion System Inlets. AIAA Paper 75-1186, Oct. 1975. 

2. Hess, J.L.: Calculation of Potential Flow about Arbitrary Three-Dimensional 

Lifting Bodies. Final Tech. Report. McDonnell Douglas Rept. No. MDC J5679-01 , 
Oct. 1972. (A somewhat condensed version is contained in Comp. Meth. in 
Appl. Mech. and Engg., Vol. 4, No. 3, Nov. 1974.) 

3. Lieblein, S.; and Stockman, N.O.: Compressibility Correction for Internal 

Flow Solutions. J. of Aircraft, Vol. 9, No. 4, Apr. 1972. 

4. Albers, J.A.: Theoretical and Experimental Internal Flow Characteristics 

of a 13.97-Centimeter-Diameter Inlet at STOL Takeoff and Approach Conditions. 
NASA TN D-7185, Mar. 1973. 

5. Mack, D.P.: Calculation of Potential Flow about Arbitrary Three-Dimensional 

Lifting Bodies. Users Manual. McDonnell Douglas Report No. MDC J5679-02, 

Oct. 1972. 

6. Halsey, N.D. ; and Hess, O.L.: A Geometry Package for Generation of Input 
Data for a Three-Dimensional Potential -Flow Program. NASA CR 2962 
(McDonnell Douglas Rept. No. MDC J7710), Mar. 1978. 

7. Halsey, N.D.: A Three-Dimensional Potential -Flow Program with a Geometry 
Package for Input Data Generation. NASA CR 145311 (McDonnell Douglas Rept. 

No. MDC J7733), Mar. 1978. 


88 



APPENDIX A. OVERLAY STRUCTURE 
(IBM 370) 


ENTRY PlAxw 
CJ V E K L A Y ALPHA 

INSERT Kct" AN» utJLNPI yPANcLCyPANcLS 
1 NsERi i KARSP ,SP ALEC* sPAL ES » LbPLN2 yCUePLN 
INjlh I AK *.LLN y Dt RV » IN t EKP y uEKV2 y i NTRP2 » NORM y AM AXA 
INSERT PLGSu»l LOSi,T LGS2,PTSI f PTSi,DERC r UERS 
INSERT SS P y SEo y P UNCH y MAXP * PLANE 3 
OVERLAY ALPHA 
INSERT jRPuT 
1 1 \‘ S E K i PL' 1 N i S 
i NSc..K i A i li'u 
OVERLAY bLI A 
INSERT RUL iP I 
OVtKLAY otU 
JlNSERI' Llt-T 
OVERLAY Alpha 
INSERT vplKPs 
inslki Rear 

lNSLkl WiiiLAK 
INScRT KOHLtVl 
INSERT fi ELDS 
OVERLAY del A 
j.NStK'1 Vr.HNLP 
UVlRl.AY lit I a 
INSERT kFKlR 
UVeKlay Alpha 
INS CH I PSWISJ: 

OVERLAY ALPHA 

INSERT API kM 
OVERLAY ALPHm 
J.NSt k« LuLSwL 
INScRl SPACth 
OVERLAY ALPHA 
i l\ S L K I p SluIN 
OVERLAY ALPHA 
INSERT lSGLOM 
OVLRlAY alpha 
JUNSlRT, LRuSLL 
OVERLAY ALPHA 
J.NSLk'1 PLPASG 
OVlKlAY ALPHA 
I. NS b R f L 0 Pi P RO 
OvLKlaY ALPHA 
INSERT ‘PiLlKi 
i N‘S Etv ( INI SL { 

i NS t K I ON Sl KN » Cf PI iSi b y S L i 3S R 

i NS 2 KT At A LL , i GS t xGSP Ll , PER S PL , P i C TR2 , ROT AT V , S L DR AW , ST 1 T L 0 

INSERT ST IT Li 

iNSfc* I ANGoLK « PLR »S iNCOS 

INS LkT lakTLS » lh LlK. » KYHUu y Re ADO y SLTRAC y TRACE » PVPIND » WRIT&W 
i NS L K « X 1 h i A. P 
A N S ER I i,ln i Ay SLS I KI 
u V jb K L A V A LPhA 
I NSc rt 1 LDrtPkS 



APPENDIX B. EXTERNAL UNITS USED IN THE INLET PROGRAM 


A total of 30 temporary storage external units are used. The unit numbers 
are: 1, 2, 3, 4, 8, .9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 
24, 25, 26, 27, 28, 29, 30, 31, 32, 33. In addition to these units, a SD-4060 
unit is also required if graphical plot option is selected. 

There are 3 modes of operation in this program: 

Mode 0: To obtain the fundamental solution, cross-section and combined 

solution (incompressible or compressible). 

Mode 1: To obtain the fundamental solution. Basic geometry data is 

input only and unit numbers 4, 22, 23, 24, 29, 30, 32, 33 should 
be saved so that data may be used for the combination solution. 

Mode 2: To obtain the combination solution only. Cross-section data and 

the combination data should be input but not the basic data that 
defined the inlet. 

Example of the external units set up is illustrated below. 
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B.l External Unit Set-Up 
(IBM 370) 


/GO. 

// 

//GO. 

// 

//GO. 

// 

//GO. 

// 

//GO. 

//GO. 

// 

//GO. 

// 

//GO. 

// 

//GO. 

// 

//GO. 

// 

//GO. 

// 

//GO. 

// 

//GO. 

// 

/ / GO • 

// 

//GO. 

// 

//GO. 

// 

//GO. 

// 

//GO. 

// 

//GO. 

// 

//GO. 

// 

//GO. 

// 

//GO. 

// 

//GO. 

// 

//GO. 

// 

//GO. 

// 

//GO. 

// 

//GO. 

// 

//GO. 

// 

//GO. 

// 

//GO. 

// 

//GO. 

// 

//GO. 

// 


TO it* GO I 

> F T G 2 P 3 C 1 

rTOirUCl 

F 1 u s i - u u i 

r T G 7 r 3 3 i 
P T u o P C 3 i 

F ( G vF G 0 l 

P T i G r U 0 1 

FT 1 IF OGi 

F I i A F o' G 1 

FTibPuGi 

t-Tlnr-GUi 

P'TlSP gw 1 

r T 1 6 P G C 1 

F T i 7 F G G i 

FTitroC i 

P T l S F GUI 

F 1 a 0 r 3 gi i 

F 7 a 1 P v. G l 

r 7 2 2 r G G i 

F T £ 2 F u o I 

FIaa-FCUa 

FTZSF.Ui 

rTZ.br 3U i 

F l A 7 r 3 0 1 

F T 2 L r G G 1 

FTtVFi-Gl 

rTbCrOOl 

F 7 2 i r OGi 

F1i>2rGCi 

FT 3 31 001 

S Oh Gc> G 


U O D S iN A i v i t. — 1« F l o i 
DCb* ( KhGFfT-VBj 
DU DSNAME=3pfG2 
DCS* ( kt:CPM=VbS 
UD U S> » A M t = b r T 0 3 
0C6 - ( ktCFM-Vo j 
DU D SN AM 2 = « F f 04 
DC b = ( R c G F iT - V B G 
3 3 GUM M Y 
O J 0 SN AMF. *tP i 'G c 
0 G 3 = l R c G F M= V b 0 
DO DS<NAKE = tFT Jy 
D L 3 — ( KL 3 F M=V 0 S 
Du DSNAME = CFi i 0 
OGb - l RE 3 I’M- V o G 
• 0 D DSN AM r = G F I ll 
DCb = tRfcGFH=VF,S 
DD U SNA Mi: = LPT la 
D Cb * ( K L C r M* V b -> 
DO OSNAMc-aFT ib 
DCb - ( RE G FiT=Vb S 
00 U SN AM f. = GF T 1 h 
D C d * l KG 3 f M — V b b 
DU D SNAMii-bFT 1 b 
DC b = l Kb G FM-V 3 3 
00 0 S NAME - OF T i 6 
DCS- ( RLClM=Vb3 
DD DSN AMD- OFT 17 
_ OCb* ( R L 3 F H =V b S 
DD _ D , SNaMc = (.F 1 it 
D Gb= ( RE C F M = V b 3 
DD DS.NAMF.-tFT L V 
G‘ 3 b= l Kh 3 H*j=V b 3 
u'J DSNmME-LPT ZU 
DGb - ( kc.CFtf=VbS 
DD DSNAMc=tF F. i 
OCB= (RCGFH-VbS 
uu DS!NAMl*GPT 2 2 
03b- ( RECrtf=VbU 


DO 


SriSl AM c - 3 r f 4 3 


DCb* l RfccFM-VBb 
DO U SNAMt = Or T 2 a 
Dob— 1 K g 3 P »'i — V 3 3 
0 0 DSN AM c — t. F T 23 
033= l kc CPH=VbS 
Ou 3 SNaML — tF I c c 

33b - ( KFCrM = Vb3 

Dp DSN aMf-uFI 2 i 

DG G - ( Rt: G F rj - V b S 
0 D G Su AM t = c F f a L 


DGb 
DD 


( KgGG V bG 


0 SN At* G — C F 7 2 V 
DCb = ( KtLFMfVbb 
D 3 3 S N /> i i i E — t F 1 3 0 
DCb - ( '<GCt"M-Vbo 
DD uSMA.H:-ur I 3 1 
OCb = t KC 3 F.T= vB G 
0 -l 3 S.N AMG — i.»F I 3 2 
Cub - ( R: Gt M=VD3 

00 OSNAMu - gF’I 30 
, t = l R c G r M - V B ! 


U.N i T i S Y S D A » S FA C L = ( TRK 
BLKS1ZF = ohh 7 tbUFNG^l ) 
UN l T = S YSO A » SF ACC- ( TRK. 
BLKSi It-bH+l ,bUf N0=1 ) 
UNI l = SYSDA tSFACcM TRK 
bLKbi ZF = 6 hh? t BUrNO=l ) 
UN i T= SYSLA } SPA Cfc= ( TRK 
b g K S i Z 0 = 6 A h 7 t B.UF NO-i ) 


Ui\I r=SYS3A t SF’ACE= ( TRK 
bLKOI Zt=6AA7 »boFNO=i ) 
UN1T=SYSDA, SPACE- ( TRK 
BL K S i Z E= i 50 3 0 » DO F NU= 1 
UN i 1 = SYSDA » SPACE = l TRK 
bLKS.lZt — o-» h 7 r b U F ii 0 — i ) 
UNi T= SYSDA, SPA Cc={TRK 
bL KS i Z fc = 1 3 03 0 , 3UF N0= 1 
UiT i T = SY SDa ,S PAC t = ( TRK 
OlKsI i l-i 3O3C , bUF NU = i 

uni t=sysda, space = i tkk 

b L K S i Z L— o h h / ,?iUFNC!=l ) 

i i\> 1 l - 3 * 3 n a .Cisrc- / 


iLNiitt-ONH 1 , bUMNU-i I 

UN i T = S Y S 0 A , SP A C E - ( IRK 
■ OLKS i ZC— 6 hhY ,BUFN 0~1 ) 
UNIT =SYSDA,SFACt:= ( TKK 
3L>\3 I L E = 6-«A 1 1 bUF.NU= 1 ) 
UNIT -SYSDA ,SPAC2= ( TRK 
BLKS I Z l.= fcAH 7 , 3UFNu=l ) 
UNIT -SYSDA , SPACE = ( TRK 

jitc-oSi / ibUFtNO— i ) 
UNi f=SYSDA,SPACF=l TRK 
b L K S i Z 0= 6 h*i 7 , a U F N 0 = i) 
UN iT« SYSDA, SPACE* ( TRK 
bGK jl ih-O'iM / ,bUPNG=i T 
UNIT -SYSOA,SPACE= (Tkk 
BLKS iZt = 6H*«7 ,bU t *ND=i T 

UNi r= sysda, space* { irk 
bLKbiZt — 6 hh 7 , BUFNlD-i ) 
UNiT-SYSOA, SPACE* l TRK 
bLKSi <_c-6W» ,;iUFNO=l ) 


- 0 3 S A - R 3 3 T 3 . C M 

D i . i I =■ T A P C i 6 , L A b b 3 



( 30 

(30 

(DO 

IOC 


(S3 

( so 

{ i»C 
( SO 
( SO 

( so 
( so 
( so 

(SO 
(3C 
( 20 
( 3C 
( 30 
(30 
(30 
(30 
(30 

i 10 
( 1C 
l ID 


5) ) , 

b) ) , 

S) ) , 
3) ) , 

5) ) , 
S) ) , 
S) ) , 
S ) ) , 
S) ) , 


S) ) , 


S) ) , 
S) ) , 
S) ) , 
S) ) , 
S) ) , 
S) ) , 
S) ) , 
S) ) , 
S) ) , 
s) ) , 
S) ), 
S) ) , 
s) ) , 
S) ) , 


. C 3M 1 . 

= k fc'T P 3* 1 


( 10. 5) ) , 

(10. 5) ), ' 

( 1 0. S ) ) , 

( it, S) ) , 

( 10 , S) ) , 
l 1 0* S ) ) , 

, OISP* ( NEW »KE 
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Unit Unit 
Variable Number 


B.2 Data Assignments to External Storage 


N0UT 

IGE0M 

NKUTTA 

NAIJ 

NRSIDE 

NT 

ND 

NSIG 

NM 

NAMIJ 

NRM 

NSETS 

NSETV 

N0FSAV 

SCSAVE 

NFLUX 

NV 

KUNIT 

IUA 

IUB 

IUNITI 

IUNIT 


Unit #1 
Unit #2 

Unit #3 

Unit #4 

Unit #8 

Unit #9 

Unit #10 

Unit #11 
Unit #12 

Unit #13 
Unit #14 
Unit #15 
Unit #16 

Unit #17 

Unit #18 
Unit #19 

Unit #20 

Unit #21 

Unit #22 

Unit #23 

Unit #24 

Unit #25 
Unit #26 

Unit #27 

Unit #28 

Unit #29 

Unit #30 
Unit #31 

Unit #32 

Unit #33 


Unit Function 
scratch units 

individual vorticity onset flows 

geometric quantities for each input panel (on-body + wake) 
intermediate (left and right) vorticity flows 
"plus" A-. matrix 

* J 

right-hand sides for "plus" matrix 
"plus" on body + off-body V.. 

* J 

scratch units to be used in C0LS0L (matrix solutions) 

source density solutions for each onset flow 
"minus" on-body + off-body V.. 

* J 

"minus" A., matrix 

* J 

right-hand side for "minus" matrix 
static (sunmed) source densities 
static velocity distributions 
scratch unit 

generated cross-section data from the input cross-section points 
volume flux data 

combined velcoities at cross-section control points 

important data are saved on this unit for restart capability 

data for streamline tracing program 
data for streamline plots 

viewing angles for the input geometry plots 

coordinates of input on-body points 
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APPENDIX C. SAMPLE INPUT AND OUTPUT FOR A SMALL PANEL NUMBER TEST CASE 


C.l Input for the Test Case 

(Card Column) 

01234567 

123*567bYl;l*3*i>b7tiV012>Oo7890l2.i*‘>t»7a9Ca23H6t>7t>90123H:>o7o90l23*567a9012 


11101*1510010 
INLET TEST CASE wilH CENTtRbGOY 1 HAKfc 1 PL. SYM. TYPE C ONE 

TfcSl lOOOllOiO FP7P 6.0 


6 I 

3 


1.0 

0. 

G. 2 

1.50 

0. 

0.2 

2.00 

0. 

0. j5 

4.00 

0. 

0.55 

1.0 

0. 

0. I 

1.50 

-0. 173205 

O.lO 

2.00 

—0. 303 i09 

0.175 

4.00 

—0.303 109 

0.175 

1.0 

0 . 

0. 1 

1.50 

-0.173205 

-G.iO 

2.00 

-G.50.3J.GV 

-0.175 

4.00 

-o.30j>ioy 

-0.1/5 

1.0 

0. 

0. 1 

1.50 

0. 

-0.2 

2.00 

0. 

-0.35 

4.00 

0. 

"0.55 

4.UOGoOO 

G.COouOi 

1.0000U021 

2.000 000 

0 • 0CK*J5l 

i . C 00 J 

0.500000 

U»000uji 

1 .^uoOOO 

0.0 

0 *00000 1 

1 .7000 * G 

0.500000 

0. G0G0G1 

2 • 4L OOOOO 

2.000033 

O.CjC^uZ 

2.460000 

4.bGQuG0 

0.000 vGZ 

xc. .nOOUOO 

6.0OUOU J 

0. OOGoOZ 

4..hOwOjG 

4.00U0u0 

-0.bocuZ5 

0.500000 1 

2.00GG00 

-v>. ct>6 0^5 

0.5v00 jU 

0.50GGG0 

-1 .C5V250 

O.bbOOoO 

0.0 

-1.472*42 

0 • 0 5 0 - v. ; v» 

0.500000 

-i.VC5<i55 

1.GV9WV 


2.00000U -2. 07 oho 1 1.199999 

4.00OuU0 — ^.GVbno 1 i.iVVVVV 
6.000C00 -2.0764oi l.lVVVVv 

4.000000 — G.bo6Gc5 -G.50O0001 

2.0GD0G0 -G.bt>6025 -G.50vK>uu 
0.500000 -i.G3V*:30 -u.oGOOuJ 
0*0 — 1.^722 4 *2 -0.t>5UUOL 

G.50CJ00 — l.Vo5255 — i.UVVvW 
2.00GJCG -2.C7b*ol -Ul^WVV 
4.000GC0 -2.07b4t>l -l.iVV^VV 
O.GuJD^O -2.67bHt>i -i .1VVVVV 
4»00u00t> -O*0ucoi»o — l.UOuOuGl 
2.00G./DO — O.GCGooO. — l.«JOOGu\> 

O.bGGuOO “0 • C JvOc v “i • ^L'OGcu « 

0,0 -O.GGOCOl -1.700000 

O.5GGQG0 — 0 • CUOGu 1 -2.2G00O0 
2 * 0 ^ 30^0 **i)» C l 1 o ^'. 1 1 **/ 0 ^ 

4.QuUuOU — O.GOotGi — 2 .hOCGG 0 
o'UC'ji "i. •**OOOJw3 

8 

0*0 G.G C. 0 VlLW GF bUDY LOURING GOWN Thl +X AXIS TOWARDS THE ORIGIN 

— 90.0 0.0 u. U View UF u JU Y lLUMi< 0 UUrtN IHc ♦Y AXIS TOwARDb l tit ORIGIN 

0*0 90*0 0. 0 VlLW OF oUUY LUUXJ.UG OUwfe The — Z AXiS IuWARuS Ilit ORIGIN 

5*0 3.0 G.C C 5 UtGRtt VlLW uF lHt t>OL>Y FROM lnt + X ♦Y +2 SlDt 

10. G 10*0 G.C* 0 10 DeGRct VlLW uF Tut tiulY FkuM Trie +X + Y **-Z SIDE 

15.0 15.0 G.u G 15 GeGRct VitW OF Iht bUUY FROM THfc +X ♦Y +2 SlOt 

30*0 30*0 0*0 G 3u GtGRtL VlcW OF FtiL LuGY FROM Tht+X +Y +Z 5 1 0 1 


-45.0 

45.0 G. 

1 45— OtbKcL 

VitW OF 

lht bub Y 

FKUH THt +X +Y *Z SiOt 

0.25 

0. 

0. 

1.0 

0. 

0. 5 

3.0 

0. 

0. 

1.0 

0. 

0. 5 

0.125 

0. 

0. 

1.0 

0. 

0. 

1.250 

0. 

0. 

1.0 

0. 

0. 

2.500 

0 . 

0. 

1.0 

0. 

0. 

3.5 00 

0 . 

0. 

1.0 

0. 

0. 

G.O 

0.0 

141.6 



1 






15.874 

0.02 

1.0 

1.0 

0.02 



0.02 

1.0 

1.0 

0.02 
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VO 


CASE NO. TEST 


PROGRAM J1HP # DOUGLAS AIRCRAFT COMPANY P AQE 

LONG BEACH DIVISION 
THURSDAY. APR L9 » 1979 

INLET TEST CASE WITH CENTERBOOY 1 WAKE l PL. SYM. TYPE C 


ECT IONS - 
POINTS 


CASE ID 

NUMBER OF LIFTING S 
NUMBER OF OFF-BCDY 
BASIC BODY INPUT CHECKOUT 

SPECIAL LAST WAKE OPTION 

PIECEWISE VORTICITY OPTION 

SYMMETRY INPUT CPTION 

DIPOLE STRENGTH CONSTANT INPUT 
FUNDAMENTAL SOLUTION OPTION — 

CROSS-SECT ION INPUT OPTION 

FLOW PASSAGE CPTION 

COMBINATION PROGRAM OPTION 


TEST 

b 

0 
I 
0 
1 
0 

1 

0 


INPUT CONTROL STATION NO. 1 

INPUT CENTERBOOY DOWNSTREAM OF INLET LIP. TYPE C 

USER SELECTS TO RUN THE FOLLOWING 

FUNDAMENTAL SOLUTION PROGRAM 
CROSS-SECT I CN PROGRAM 
FLOW PASSAGE PROGRAM 

COMB I NAT I ON PROGRAM ( COMPRESSIBLE ) 
COMPONENTS OF THE UNIFORM ONSET FLOWS 
III l. OOOOOO. 0.0 , 0.0 

I 2 ) 0.0 . 0.0 . 1.000000 

I 3 ) 0.0 , 1.000000, 0.0 

INPUT SCALE FACTOR VALUE 6.000000 


ONE 


* — 


C.2 Output for the Test Case 



PROGRAM JIHPt DOUGLAS AIRCPAFT COMPANY PAGE 2 

LONG BEACH DIVISION 
CASE NO. TEST THURSDAY f APR 19, 1979 

INLET TEST CASE WITH CENTERBODY 1 WAKE 1 PL. SYM. TYPE C ONE 


N 

M 

X 

X 

X 

X 

NX 

XO 

D 

TYPE OF 



Y 

Y 

Y 

Y 

NY 

YO 

T 

ELEMENT 


l Z JL l NZ _ZO __ A 


X 

1 

6.000000 

0.0 

0.0 

5.000000 

0.0 

1.200000 

9.000000 

-1.039229 

0.600000 

6.000000 

0.0 

0.0 

-0.327327 

-0.472456 

0.818317 

7.999998 

-0.346410 

0.600000 

2.3842E-07 
3 • 2 31 IE + 00 
1 .9049E+00 

NL IF 


2 

9.000000 

0.0 

1.200000 

12.000000 

0.0 

2.099999 

12. 000000 
-1.818653 
1.049999 

9.000000 

-1.039229 

0.600000 

-0.251459 
-0. 48 3934 
0.838198 

10.636363 

-0.732185 

1.268181 

2 • 0489F— 07 
3.511 4E + 00 
5.1143F+00 



3 

12.000000 

0.0 

2.099999 

24.000000 

0.0 

2.099999 

24. 000000 
-1.318653 
1.049999 

12.000000 

-1.818653 

1.049999 

0.0 

-0.500000 

0.866025 

18.000000 

-0.909327 

1.574998 

5 • 3644E-07 
1 .2 182E+01 
2.5200E+01 


2 

1 

6.000000 

0.0 

0.0 

9.000000 

-1.039229 

0.600000 

9.000000 

-1.039229 

-0.600000 

6.000000 

0.0 

0.0 

-0.327326 
-0.9449 11 
0.0 

7.999998 

-0.692819 

0.000000 

l .192*16-07 
3.231 IE + 00 
1 • 9049E +00 

NLIF 


2 

9.000000 

-1.039229 

0.600000 

12.000000 

-1.818653 

1.049999 

12.000000 
-1.818653 
-1. 049999 

9.000000 

-1.039229 

-0.600000 

-0.251460 

-0.967868 

0.0 

10.636362 

-1.464368 

-0.000000 

3 . 5 763F-07 
3 .5114E+00 
5.1 143E+00 



3 

12.000000 

-1.318653 

1.049999 

24.000000 
-1.818653 
l. 049999 

24. 000000 
-1.818653 
-1. 049999 

12.000000 

-1.818653 

-1.049999 

0.0 

-1.000000 

0.0 

18.000000 
-1.8 18653 
0.0 

0.0 

1 .2182E+01 
2 • 5200E+01 


3 

1 

6.000000 

0.0 

0.0 

9.000000 
- 1.039229 
-0.600000 

9.000000 

0.0 

-1.200000 

6.000000 

0.0 

0.0 

-0.32 73 2 7 
-0.472456 
-0.818317 

7.999998 

-0.346410 

-0.600000 

2.3842E-07 
3.2311F +00 
1.9049E+00 

NLIF 


2 

9.000000 

-1.039229 

-0.600000 

12.000000 

-1.818653 

-1.049999 

12.000000 

o.o 

-2. 099999 

9.000000 

0.0 

-1.200000 

-0.251459 

-0.483934 

-0.838198 

10.636363 

-0.732105 

-1.268180 

0.0 

3.5U4E + 00 
5.1143E+00 



3 

12.000000 

-1.818653 

-1.049999 

24.000000 
-1.818653 
- 1.049999 

24.000000 

0.0 

-2. 099999 

12.000000 

0.0 

-2.099999 

0.0 

-0.500000 

-0.866025 

18.000000 
-0.909326 
- 1.574998 

2 .9802E-07 
1.2182F+01 
2.52006+01 




************ 


************ 

************ 


1 

1 

24.000000 

0.000006 

6.000000 

12.000000 
0. 000006 
6.000000 

12. 000000 
-5. 196150 
3.000000 

24.000000 
-5. 196150 
3.000000 

0.0 

0.500000 
-0. 866026 

18.000000 

-2.598071 

4.500000 

2.4000F+01 
1.3416E+0I 
7.2 OOOE + 01 

LIFT 


2 

12.000000 

0.000006 

6.000000 

3.000000 
0.000006 
7. 199999 

-l 1:852888 

3.599999 

3.000000 

-Q. 11 4708 
0.496699 
-0.860309 

7.500000 

-2.857880 

4.949999 

l .8120E+01 
1.1225E+01 
5.9926E+01 



cn 



VO 

CD 

PROGRAM JIHP. DOUGLAS AIRCRAFT COMPANY PAGE 3 

LONG BEACH DIVISION 
CASE NO. TEST THURSDAY, APR 19, 1979 

INLET TEST CASE WITH CENTERBODY 1 WAKE 1 PL. SYM. TYPE C ONE 


N 

H 

X 

Y 

Z 

X 

Y 

Z 

X 

Y 

Z 

X 

Y 

Z 

1 

££21 

XO 

YO 

zo 

0 

T 

A 

TYPE OF 
ELEMENT 

1 

3 

3.000000 

0.000006 

7.199999 

8l §00006 

10.199999 

-si 333454 
5.099999 

3.000000 

-6.235382 

3.599999 

-0.654653 

0.377964 

-0.654654 

1 . 50C000 
-3*767204 
6.524998 

7*93738+00 
9*67948+00 
3.691 IE+01 

LIFT 


4 

8*So0006 

10.199999 

oi §88886 

13.199999 

-i?: ISMSS 

6.599997 

0.0 

-8.833454 

5.099999 

-0.654654 
-0. 377965 
0.654654 

1.500000 

-5.066241 

8.774994 

7.2373EV-00 
1 «2445E*0l 
4.9639E»-01 



5 

3.000000 

0.000006 

13.199999 

12.000000 

0.000012 

14.399998 

12.000000 
-12.470764 
7. 199993 

3.000000 

-11.431532 

6.599997 

-0. 114707 
-0.496700 
0.860309 

7.500000 

-5.975566 

10.349991 

1.81205+01 
1 .651 96+01 
1 .25306+02 



6 

12.000000 

0.000012 

14.399998 

24.000000 

0.000012 

14.399998 

24. 000000 
-12.470764 
7.199593 

12.000000 

-12.470764 

7.199993 

0.0 

-0.500000 

0.866026 

18.000000 

-6.235375 

10.799988 

2. 40006+01 
1.87456+01 
1 .72806+02 



7 

2 A. 000000 
0.000012 
14.39999B 

36.000000 

0.000012 

14.399998 

36. 000000 
-12.470764 
7. 199993 

24.000000 

-12.470764 

7.199993 

0.0 

-0.500000 

0.866026 

30.000000 

-6.235375 

10.799988 

2.40006+01 
l .87456+01 
1 .7230E +02 

WAKE 

2 

1 

24.000000 

-5.196150 

3.000000 

12.000000 

-5.196150 

3.000000 

12.000000 

-5.196150 

-3.000000 

24.000000 
-5.196150 
-3 .000000 

0.0 

1.000000 

0.0 

18.000000 

-5.196150 

0.0 

2.40006+01 
1 .34166+01 
7.2000F+01 

LIFT 


2 

12.0COOOO 
-5. 196150 
3.000000 

3.000000 

-6.235382 

3.599999 

3.000000 

-6.235382 

-3.599999 

12.000000 
-5.196150 
-3 .000000 

-0. 114708 
0.993399 
0.0 

7.500000 

-5.715766 

0.0 

1.81206+01 
1.1225E+01 
5.992 66+01 



3 

3.000000 

-6.235382 

3.599999 

0.0 

-8. 833454 
5.099999 

0.0 

-8.833454 

-5.099999 

3.000000 

-6.235382 

-3.599999 

-0.654653 

0.755929 

0.0 

1.500000 

-7.534418 

0.0 

7. 93726+00 
9 . 67946+00 
3.691 IF +01 



4 

0.0 

-8.833454 

5.099999 

3.000000 

-11.431532 

6.599997 

3.000000 
-1 1.431532 
-6. 599997 

0.0 

-8.833454 

-5.099999 

-0.654654 

-0.755929 

0.0 

1.500000 

-10.132492 

0.0 

7.93736+00 
1 • 2445E +0 1 
4.9639E+01 



5 

3.000000 

-11.431532 

0.599997 

12. 000000 
-12.470764 
7. 1 99993 

12.000000 

-12.470764 

-7.199993 

3.000000 

-11.431532 

-6.599997 

-0. 114708 
-0.993399 
0.0 

7.500000 

-11.951141 

0.0 

1 • 81206+01 
1 .65196+01 
l .25306+02 



6 

12.000300 

-12.470764 

7.199993 

24.000000 

-12.470764 

7.199993 

24. 000000 
-12.470764 
-7. 199993 

12.000000 
-12 .470764 
—7 • 1° 999 3 

0.0 

-1.000000 

0.0 

18.000000 
-12.470 764 
0.0 

2.4000E+01 
1 .87456+01 
1 .72806+02 



7 

24.000000 

-12.470764 

7.199993 

36.000000 

-12.470764 

7.199993 

36.000000 

-12.4707*4 

-7.19999J? 

24.000000 
-12.470764 
-7. 199993 

0.0 

-1.000000 

0.0 

30.000000 

-12.470764 

0.0 

2.40006+01 
1 .8745E+01 
X .72806+02 

WAKE 



PROGRAM J1HP, DOUGLAS AIRCRAFT COMPANY PAGE 4 

LONG BEACH DIVISION 
CASE NO. TEST THURSDAY, APR 19, 1979 

INLET TEST CASE WITH CENTERBODY l WAKE l PL. SYM. TYPE C ONE 

* 


M 

X 

X 

X 

X 

NX 

XO 

0 


Y 

Y 

Y 

Y 

NY 

YO 

T 


l 

Z 

Z 

Z 

NZ 

ZO 

A 

l 

24.000000 

12.000000 

12.000000 

24.000000 

0.0 

18.000000 

2.4000F+01 


-5.196150 

-5.196150 

0.0 

0.0 

0. 500000 

-2.598075 

l .3416F+01 


-3.000000 

-3.000000 

-6.000000 

-6.000000 

0.866025 

-4.500000 

7.2000F+01 

2 

12.000000 

3.000000 

3. 000000 

12.000000 

-0.114708 

7.500000 

1 • 3 12 0E + 01 


— 5.196150 

-6.235382 

0.0 

0.0 

0.496700 

-2.357382 

l . 1 2? 5 P + 01 


-3.000J00 

-3.599999 

-7.199999 

-6.000000 

0.860309 

-4.949999 

5 • 992 6E +01 

3 

3.000000 

0.0 

0. 0 

3.000000 

-0.654654 

1.500000 

7.93 73E + 00 


-6.235332 

-8.833454 

-0.000006 

0.0 

0. 377964 

-3.767209 

9.6793E+00 


-3.599999 

-5.099999 

-10. 199999 

-7.199999 

0.654654 

-6.524998 

3 • 691 IE + 01 

4 

0.0 

3.000000 

3. OOOOOO 

0.0 

-0.654654 

1.500900 

7.9373F +00 


— 3*833454 

-l 1.431532 

-0. 000006 

-0.000006 

-0.377965 

-5.066248 

1 .2 445E +01 


-5.099999 

-6.599997 

-13.199999 

-10.199999 

-0.654654 

-8.774994 

4.9639E+01 

5 

3.000000 

12.000000 

12.000000 

3.000000 

-0.114708 

7.500000 

1 • 8 120F+01 


- 11.431532 

-12.470764 

-0. 000006 

-0.000006 

-0.496700 

-5.975574 

1 .6519E+01 


-6.599997 

-7.199993 

-14.399998 

-13.199999 

-0.860309 

-10.349991 

1 .2530E + 02 

6 

12.000000 

24.000000 

24. OOOOOO 

12.000000 

0.0 

18.000000 

2.4000E +01 


— 12.4707 64 

-12.470764 

-0. 000006 

-0.000006 

-0.500001 

-6.235385 

1 .8745E+01 


- 7.199993 

-7.199993 

-14.399998 

-14.399998 

-0.866026 

-10.799988 

1 .7280E+02 

7 

24.000000 

36.000000 

36. OOOOOO 

24.000000 

0.0 

30.000000 

2*4000^+01 


-12.470764 

-12.470764 

-0.000006 

-0.000006 

-0.500001 

-6.235385 

1 .8745F+01 


-7.199993 

-7.199993 

-14.399998 

-14.399998 

-0.866026 

-10.795988 

1.7280E+02 


************ ************ ************ 


TYPE OF 


ELEMENT 


LIFT 


WAKE 


to 
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CASE NO. TEST 


PROGRAM J1HP, OOUGLAS AIRCRAFT COMPANY PAGE 

LONG • BEACH DIVISION 
THURSDAY* APR 19. 1979 

INLET TEST CASE WITH CENTERBODY 1 WAKE l PL. SYM. TYPE C 

* 


5 


ONE 


INPUT 

SECTION NO. 
1 
2 


TABLE OF INPUT INFORMATION 


SECTION 

TOTAL NO. OF ELEMENTS 

EXTRA 

STRIP 

SOURCE ELEMENTS 

WAKE ELEMENTS IN 

TYPE 

IN EACH SECTION 

STRIPS 

NO. 

IN THE STRIP 

THE STRIP 

0 

9 


1 

3 

0 




2 

3 

0 




3 

3 

0 

1 

21 


4 

6 

1 




5 

6 

1 




6 

6 

l 


PROGRAM J1HP. DOUGLAS AIRCRAFT COMPANY PAGE 8 

LONG BEACH DIVISION 
CASE NO. TEST THURSDAY, APR 19, 1979 

INLET TEST CASE WITH CENTERBODY 1 WAKE 1 PL. SYM. TYPE C ONE 

BODY SECTION NO. = 2 TYPE = l TOTAL NO. OF POINTS = 21 NO. OF STRIPS * 

LIFTING SECTION NO. I NO. OF SOURCE ELEMENTS 6 

NO. OF WAKE ELEMENTS 1 TOTAL NO. OF ELEMENTS PER STRIP 7 

TOTAL NO. OF CONTROL POINTS ( INCL. OFF BODY POINTS ) = 27 

TOTAL NO. OF ELEMENTS IN THE LIFTING SECTION = 21 

NO. OF FAR ELEMENTS = 36 NO. OF INTERMEDIATE ELEMENTS = 16B NO. OF NEAR ELEMENTS * 

THE 27 X 27 MATRIX WITH 5 RIGHT SIDES WAS SOLVED DIRECTLY IN 0 .0 1 1MINUTES. 

THE 27 X 27 MATRIX WITH 1 RIGHT SIDES WAS SOLVED DIRECTLY IN 0 .0 10M INUTE $. 


1416 



PROGRAM J 1HP f DOUGLAS AIRCRAFT COMPANY PAGE II 

LONG BEACH DIVISION 
CASE NO. TEST THURSDAY, APR 19, 1979 

INLET TEST CASE WITH CENTERBODY I WAKE l PL. SYM. TYPE C ONE 

FUNDAMENTAL SOLUTION OUTPUT 


I 

STATVX 

STATVY 

STATVZ 

UVXl 

UVY1 

UVZ1 

UVX2 

UVY2 

UVZ2 

UVX3 

UVY3 

UVZ3 

1 

16.91554 

-2.92990 

5.07462 

0.954 74 

-0. 16537 

0.28642 

0.21455 

0.32402 

0.27289 

-0. 12387 

0.64703 

0.32402 

2 

22.12704 

-2.87444 

4.97854 

1.22937 

-0.15970 

0.2 7661 

0.14612 

0.28300 

0.20723 

-0.08436 

0.53402 

0.28301 

3 

22.35988 

-0.00003 

-0.00001 

1.23063 

-0.00000 

-0.00000 

0.02573 

0.35573 

0.20538 

-0.01486 

0.61614 

0.35573 

4 

16.91557 

-5.35971 

-0. 00010 

0.95474 

-0.33073. 

-0.00000 

-0.00000 

-0.00000 

0.83410 

-0.24774 

0.08582 

-0.00000 

5 

22.12706 

-5.74877 

-0.00008 

1.22937 

-0.31940 

0.00000 

0.03000 

-0.00000 

0.69741 

-0. 16873 

0.04384 

0.00000 

6 

22.35986 

0. 00000 

-0.00005 

1.23063 

0. 00000 ' 

0.00000 

-0.00000 

-0.00000 

0. 82152 

-0. 02971 

0. 00000 

-C. 00000 

7 

16.91560 

-2.92991 

-5. 07468 

0.95474 

-0. 16536 

-0.28642 

-0.21454 

-0.32401 

0.27289 

-0. 12 367 

0.64703 

-0.32402 

8 

22.12708 

-2.87436 

-4. 97859 

L. 22937 

-0.15970 

-0.27661 

-0.14612 

-0.23300 

0.20723 

-0. 08436 

0.53402 

-0.28301 

9 

22.35986 

0.00001 

-0. 00000 

1.23063 

-0. 00000 

0.00000 

-0.02573 

-0.35573 

0.20538 

-0. 01466 

0.61614 

-0. 35573 

10 

22.45866 

-0.00001 

0. 00000 

1.19751 

-0. 00000 

0.00000 

0.02409 

0.33969 

0.19613 

-0.01 391 

0.58837 

0.33969 

11 

20.20905 

1.16665 

-2. 02096 

1.13273 

0. 06540 

-0.1 1327 

-0. 14185 

0.32269 

0.20523 

0. 08190 

0.57784 

0.32269 

12 

8.04817 

3.48491 

-6. 03613 

0.33951 

0. 14701 

-0.25463 

-1 . 16699 

-0.029H 

1.15019 

0.67376 

1.11658 

-0.02910 

13 

-1.35973 

0.58878 

-1.01980 

0. 92633 

-0.40111 

0.69475 

1.32015 

-0.01053 

1.31408 

-0.76219 

1.30192 

-0.01 053 

14 

-0.35757 

0.02065 

-0. 03576 

1. 19592 

-0.06905 

0.11959 

0.34165 

0.61203 

0.39391 

-0. 19725 

1.10562 

0.61203 

15 

-2.29943 . 

-0.00002 

0.00011 

1.07329 

-0. 00000 

0.00000 

-0. 10418 

0.63284 

0.36537 

0. 10633 

1.0961 1 

0.63284 

16 

22.45863 

-0.00002 

-0. 00003 

1.19751 

-0. 00000 

0.00000 

0.00000 

0.00000 

0.78449 

-0.02782 

0.00000 

-0.00000 

17 

20.20906 

2.33352 

-0. 00002 

1. 13274 

0. 13080 

0.00000 

0.0 

0.00000 

0.76415 

0.16380 

0.01892 

-0. 00000 

18 

8.04815 

6.96 99 0 

-0.00003 

0.33951 

0.29402 

0.00000 

-0.00000 

-0.00000 

l. 05978 

1.34752 

1.16699 

-0. 00000 

19 

-1.359 75 

1.1 775 7 

-0. 00000 

0.92633 

-0.00222 

0.00000 

0.00000 

-0.00000 

1.29584 

- 1. 52437 

1.32015 

-0. 00000 

20 

-0.35756 

0.04129 

- 0. 00001 

1.19592 

-0.13809 

0.00000 

0.0 

0.00000 

1.45398 

-0.39450 

0.04556 

0.00000 

21 

-2.29944 

-0.00014 

-0.00003 

1.07329 

-0. 00000 

0.00000 

-0.00000 

0.00000 

1.46148 

0.21267 

0.00000 

0. 00000 

22 

22.45863 

0.00000 

-0. 00001 

1. L9751 

-0. 00000 

-0.00000 

-0 .02409 

-0.33969 

0. 19612 

-0. 01391 

0.58837 

-0.33969 

23 

20.20901 

1.16671 

2. 02094 

1. 132 73 

0. 06540 

0.1 1327 

0. 14185 

-0.32269 

0.20522 

0.08190 

0.57784 

-0.32269 

24 

8.04811 

3.48491 

6.03609 

0.33950 

0. 14701 

0.25463 

1 . 16699 

0.02911 

1. 1 5019 

0.67376 

1.11658 

0. 0291 L 

25 

-1.359 74 

0.58879 

1. 01979 

0.92633 

-0.40111 

-0.69675 

-1.32015 

0.01053 

1.31407 

-0. 76219 

1.30192 

0.01 053. 

26 

-0.35757 

0.02064 

0. 03577 

1. 19592 

-0. 06905 

-0.1 1959 

-0.34164 

-0.61203 

0.39891 

-0.19725 

1.10562 

-0.61203 

27 

-2.29945. 

-0.00004 

-0.00008 

1.07330 

-0. 00000 

-0.00000 

0 . 18418 

-0.63284 

0.36537 

0. 10634 

1.09611 

-0.63284 


VO 

VO 
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PROGRAM JIHPi 


CASE NO. TEST 


DOUGLAS AIRCRAFT COMPANY PAGE 12 

LONG REACH DIVISION 
THURSDAY* APR 19, 1979 


INLET TEST CASE WITH CENTERBQOY 1 WAKE 1 PL. $YM. TYPE C 
CROSS SECTION AND CONTROL STATION PROGRAM 


CROSS SECTION POINT 
YCSP 


0# 150000E + 01 
0.1 8000QE+02 


0,100000 E*Oi 
0.1000006*01 


NORMAL VECTOR 
UJ 


FORMATION OF CROSS SECTION ELEMENTS 


CROSS SECTION PO TNT 0.1 
ELEMENTS TO BE FORMED PER ! 
NORMAL VECTOR COMPONENTS 
UNIT NORMAL VECTOR 
All A 12 A 13 
A21 A22 A23 

A3 1 A32 A 33 


;^ooce+oi 5 


0. iOOOOOE+Oi 

o.ioooooe+oi 

0.0 

0.0 

0. IOOOOOE + OI 


CROSS SECTION NO. 


0.0 

0.0 

0. 1 OOOOQE+O 1 
0.0 
0.0 


LOOOOO E+0 1 


in in 



PROGRAM J1HP 


CASE NO. TEST 

INLET TEST CASE WITH 
CROSS SECTI 



X 

X 

X 


Y 

Y 

Y 


l 

Z 

1 

1 

1.50000 

1.50000 

1.50000 


0.0 

0.00000 

-1.50688 


0.0 

1.74000 

0.87000 

2 

1.50000 

1.50000 

1. 50000 


0. 00000 

0.00000 

-3.01377 


l. 74000 

3.48000 

1. 74000 

3 

1.50000 

1.50000 

1.50000 


0.00000 

0.00000 

-4. 52065 


3.43000 

5.22000 

2. 61000 

4 

1.50000 

1.50000 

-him 


0. 00000 

0.00000 


5.22000 

6.96000 

3.43000 

5 

1.50000 

1.50000 

1.50000 


0. 00000 

0.00001 

-7. 53442 


6. 96000 

8.70000 

4.35000 


1 

1.50000 

0.0 

0.0 

1.50000 

-1.50688 

0.87000 

1. 50000 
-1 . 50638 
-0.87000 

2 

1.50000 

-1.50683 

0.87000 

1.50000 

-3.01377 

1.74000 

1.50000 
-3. 01377 
-1. 74000 

3 

1.50000 

-3.01377 

1.74000 

1.50000 

-4.52065 

2.61000 

1. 50000 
-4. 52C65 
-2.61000 

4 

1.50000 
-4. 52065 
2.6L090 

1.50000 

-6.02753 

3.48000 

1.50000 
-6. 02 753 
-3.48000 

5 

1. 50000 
-6. 02753 
3.48000 

1.50000 

-7.53442 

4.35000 

1. 50000 
-7. 53442 
-4.35000 


1.53000 

0.0 

0.0 

1.50000 

-1.50688 

-0.87000 

1.50000 
-0.00000 
-1. 74000 



X 

RIX 

R2X 

xc 

XI 

TYPE 

Y 


R2Y 

*5 

ETA 



li 


1.50000 

0.0 

0.0 

1.50000 

-0.37672 

CROS 

0.0 

-1.50688 

-0.00000 

-0.37672 

0.65250 


0.0 • 

0.87000 

-1.74000 

0.65250 

1.31095 


1.50000 

0.0 

0.0 

1.50000 

-1.13016 

CROS 

-1.50688 

-3 .01377 

-1.50688 

-1.13016 

1.95750 


0.87000 

0.0 

-2.61000 

1.95750 

3.93296 


1.50000 

0.0 

0.0 

i.soodo 

-1.89360 

CROS 

-3.01377 

-4.52065 

-3.01377 

-1.88360 

3.26250 


1.74000 

-0.87000 

-3.48000 

3.26250 

6.55494 


1.50000 

0.0 

0.0 

1.50000 

-2.63704 

CROS 

-4.52065 

-6.02754 

-4.52065 

-2.63704 

4.56750 


2.61000 

-1.740 0 0 

-4.35000 

4. 56750 

9.17691 


1.50000 

0.0 

0.0 

1.50000 

-3.39048 

CROS 

-6.02753 

-7.53442 

-6.02754 

-3.39048 

5.97250 


3.48000 

-2.61000 
— * 

-5.22000 

5.87250 

11. 79890 


1.50000 

0.0 

0.0 

1.50000 

-0.75344 

CROS 

0.0 

-1 .50688 

1.50688 

-0.75344 

0.0 


0.0 

-0.87000 

-0.87000 

0.0 

1.31099 


1.50000 

0.0 

0.0 

l .50000 

-2.26032 

CROS 

-1.50688 

-1.53688 

1 .5068 8 

-2 .26032 

-0.00000 


-0.37000 

-2.61000 

-2.61000 

-0.00000 

3.93296 


1.50000 

0.0 

0.0 

1.50000 

-3.76721 

CROS 

-3.01377 

-l .50608 

l .50688 

-3 .76721 

0.0 


-1.74000 

-4.35000 

-4.35000 

0.0 

6.55494 


1.50000 

0.0 

0.0 

1.50000 

-5.27409 

CROS 

-4.52065 

-1 .50688 

1.50688 

-5.27409 

0.0 


-2.61000 

-6.09000 

-6.09000 

0.0 

9.1 769 1 


1.50000 

0.0 

0.0 

1.50000 

-6.78097 

CROS 

-6.02753 

-1 .50688 

1.50688 

-6.78097 

0.0 


-3.48000 

-7.83000 

-7.83000 

0.0 

11.79890 


50000 

0.0 

0.0 

1.5000 

-0.37672 

CROS 

U.O 

-0 .00000 

l. 50688 

-0.37672 

-0.65250 


0.0 

-1.74000 

0.87000 

-0.65250 

1.31099 




o 

ro 


CASE NO. TEST 


PROGRAM JlHPt OOUGLAS AIRCRAFT COMPANY PAGE 

LONG BEACH DIVISION 
THURSDAY, APR 19, 1979 

INLET TEST CASE WITH CENTERBODY 1 WAKE l PL. SYM. TYPE C 
CROSS SECTION AND CONTRCL STATION PROGRAM 


L6 


ONE 


3 2 

3 


X 

X 

X 

X 

R IX 

R2X 

XC 

XI 

Y 

Y 

Y 

Y 

RiY 

R2Y 

YC 

ETA 

Z 

Z 

Z 

Z 

R1Z 

R2Z 

ZC 

AREA 

1.50000. 

1.50000 

1.50000 

1.50000 

0.0 

0.0 

1.50000 

-1.13016 

-1. 50688 

-3.01377 

-0.00000 

-0.00000 

1.50688 

3.01377 

-1.13016 

-1.95750 

-1.95750 

-0.87000 

-1.74000 

-3.48000 

-1.74000 

-2.61000 

0.0 

3.93296 

1.50000 

1.50000 

1.50000 

1.50000 

0.0 

0.0 

1.50000 

-1.88360 

-3.01377 

-4.52065 

-0. 00000 

-0.00000 

3.01376 

4.52065 

-1.88360 

-3.26250 

-1.74000 

-2.61000 

-5 • 22 COO 

-3.48000 

-3 .48000 

-0.87000 

-3.26250 

6.55494 

1.50000 

1.50000 

1. 50000 

1.50000 

0.0 

0.0 

1.50000 

-2.63705 

-2.63705 

-4.52065 

-6.02753 

-0. 00000 

-0.00000 

4.52065 

6.02753 

-4.56750 

-2.61000 

-3.48000 

-6. 96000 

-5.22000 

-4.35000 

-1.74000 

-4. 56750 

9.17691 

1.50000 

1.50000 

1.50000 

1.50000 

0.0 

0.0 

1.50000 

-3.39049 

-6.02753 

-7.53442 

-0. 00000 

-0.00000 

6.02753 

7.53442 

-3.39049 

-5.87249 

-3.48000 

-4.35000 

-8. 70000 

-6.96000 

-5.22000 

-2.61000 

-5.87249 

11.79889 


FINISHED COMPUTATION FOR CROSS SECTION 
TOTAL NUMBER OF OFF-BODY POINTS GENERATED = 


15 


CROSS SECTION POINT 


0.180000E+02 0.0 


ELEMENTS TO BE FORMED PER STRIP 
NORMAL VECTOR COMPONENTS 0.100000E+01 

UNIT NORMAL VECTOR G.100000E+01 

All A 12 A 13 0.0 

A21 A22 A23 0.0 

A3 1 A 32 A33 O.IOOOOOE+Oi 


FORMATION OF CRCSS SECTION ELEMENTS 

CROSS SECTION NO. 2 

0.0 

0.0 
0.0 

0. 1 OOOOOE+O 1 
0.0 
0.0 


0.0 

0.0 

0.0 

0.1 OOOOOE+O 1 

0.0 


TYPE 

CROS 

CROS 

CROS 

CROS 



PROGRAM J1HP, DOUGLAS AIRCRAFT COMPANY PAGE 17 

LONG BEACH DIVISION 
CASE NO. TEST THURSDAY, APR 19, 1979 

INLET TEST CASE WITH CENTERBODY 1 WAKE 1 PC. SYM. TYPE C ONE 

CROSS SECTION AND CONTRCL STATION PROGRAM 


J 

I 

X 

Y 

Z 

X 

Y 

Z 

X 

Y 

Z 

X 

Y 

Z 

R1X 

R1Y 

R1Z 

R2X 

P2Y 

R2Z 

xc 

YC 

zc 

XI 

FTA 

AREA 

TYPE 

1 

1 

18,00000 

0.0 

2.10000 

18.00000 

0.00000 

2.88000 

i 5:282?$ 

1.44000 

1.05000 

*49415 

- 0.66000 

0.0 

- 1.81865 

- 1.83000 

18.00000 

- 1.07820 

1.86750 

- 1.07820 
1.86750 
1 .68199 

CROS 


2 

18.00000 

0.00000 

2.88000 

18.00000 

0.00000 

3.66000 

18.00000 
- 3. 16965 
l. 83000 

18.00000 

- 2.49415 

1.44000 

0.0 

- 3.16965 

- 1.05000 

0.0 

- 2.49415 

- 2.22000 

ie. ooooo 

- 1.41595 

2.45250 

- 1.41595 

2.45250 

2.20888 

CROS 


3 

18.00000 

0.00000 

3.66000 

18.00000 

0.00000 

4.44000 

18.00000 

- 3.84515 

2.22000 

18.00000 

- 3.16965 

1.83000 

0.0 

- 3.84515 

- 1.44000 

0.0 

- 3.16965 

- 2.61000 

18.00000 
- 1. 75370 
3.03750 

- 1.75370 

3.03750 

2.73577 

CROS 


4 

18.03000 

0.00000 

4.44000 

18.00000 

0.00000 

5.22000 

18. 00000 
- 4. 52065 
2.61000 

18.00000 

- 3.84515 

2.22000 

0.0 

- 4.52065 

- 1.83000 

0.0 

- 3.84516 

- 3.00000 

18.00000 

- 2.09145 

3.62250 

- 2.09145 

3.62250 

3.26266 

CROS 


5 

18.00000 

0.00000 

5.22000 

18.00000 

0.00001 

6.00000 

18.00000 

- 5.19615 

3.00000 

1 8.00000 
- 4.52065 
2.61000 

0 .0 

- 5. 19615 
- 2.22000 

0 .0 

- 4.52066 

- 3.39000 

18.00000 

- 2.42920 

4.20750 

- 2.42920 

4.20750 

3.78955 

CROS 

2 

1 

±!.*ra 

1.05000 

1.44000 

- 1.44 000 

18.00000 

- 1.81865 

- 1.05000 

0.0 

- 0.67550 

- 2.49000 

0 *.67550 
- 2.49000 

if :?888 
0.0 

- 2.1 5640 
0.0 

1.68199 

CROS 


2 

18.00000 

- 2.49415 

1.44000 

18.00000 

- 3.16965 

1.83000 

18 . 00000 
- 3.16965 
- 1.83000 

18.00000 

- 2.49415 

- 1.44000 

0.0 

- 0.67550 

- 3.27000 

0.0 

0.67550 

- 3.27000 

18.00000 

- 2.83190 

0.0 

- 2.83190 

0.0 

2 . 20B88 

CROS 


3 

18.00000 

- 3.16965 

1.83000 

18.00000 

- 3.84515 

2.22000 

13. 00000 
- 3.84515 
- 2.22000 

18.00000 

- 3.16965 

- 1.83000 

0.0 

- 0.67550 

- 4.05000 

0.0 

0.67550 

- 4.05000 

-!*.50743 

0.0 

- 3.50740 

0.0 

2.73577 

CROS 


4 

18.00000 

- 3.84515 

2.22000 

18.00000 

- 4.52065 

2.61000 

18.00000 
- 4 . 52 065 
- 2.61000 

18.00000 

- 3.84515 

- 2.22000 

0.0 

- 0.67550 

- 4.83000 

0.0 

0.67550 

- 4 . Q3000 

18.00000 

- 4.18290 

0.0 

- 4.18290 

0.0 

3.26266 

CROS 


5 

18.00000 

- 4.52065 

2.61000 

18.00000 

- 5.19615 

3.00000 

18 . 00000 
- 5.19615 
- 3.00000 

18.00000 

- 4.52065 

- 2.61000 

0 .0 

- 0.67550 

- 5.61000 

0.0 

0.67550 

- 5.61000 

18.00000 

- 4.85840 

0.0 

- 4.85840 

0.0 

3.78955 

CROS 

3 


18.00000 

- 1.81865 

- 1.05000 

18.00000 

- 2.49415 

- 1.44000 

18.00000 

0.0 

- 2. 88000 

; 00000 
U.O 

- 2.10000 

0 .0 

1.81865 

- 1.83000 

0.0 

2.49415 

- 0.66000 

18 . 0000 . 

- 1.07820 

- 1.86750 

- 1.07820 

- 1.86750 

1.68199 

CROS 


o 

to 



o 

-p* 


J I 


3 

4 

5 
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LONG BEACH DIVISION 
CASE NO « TEST THURSDAY, APR 19, 1979 

INLET TEST CASE WITH CENTERBODY 1 WAKE 1 PL. SYM. TYPE C ONE 

CROSS SECTION ANO CONTRCL STATION PROGRAM 


X 

X 

X 

X 

Y 

Y 

Y 

Y 

Z 

Z 

7. 

z 

18. 00000 

18.00000 

18.00000 

18.00000 

-2.49415 

-3.16965 

0.0 

0.0 

-1.44000 

-1.83000 

-3.66000 

-2.88000 

18.00000 

18.00000 

18.00000 

18.00000 

-3. 16965 

-3.84515 

0.0 

0.0 

-1.83000 

-2.22000 

-4.44000 

-3.66000 

18.00000 

18.00000 

18.00000 

18.00000 

-3.84515 

-4.52065 

0.0 

0.0 

-2.22000 

-2.61000 

-5.22000 

-4.44000 

18.00000 

18.00000 

18.00000 

18.00000 

-4.52065 

-5.19615 

0.0 

0.0 

-2.61000 

-3.00000 

-6. 00000 

-5.22000 


R1X 

R2X 

xc 

XI 

TYPE 

R1Y 

R2Y 

YC 

FTA 


RIZ 

R2Z 

ZC 

AREA 

— 

0.0 

0.0 

18.00000 

-1.41595 

CROS 

2.49415 

3.16965 

-1.41595 

-2.45250 


-2.22000 

-1.05000 

-2.45250 

2.20888 


0.0 

0.0 

18.00000 

-1.75370 

CROS 

3.16965 

3.84515 

-1.75370 

-3.03750 


-2.61000 

-1.44000 

-3.03750 

2.73577 


0.0 

0.0 

18.00000 

— 2 . 09t 45 

CROS 

3.84515 

4.52065 

-2.09145 

-3.62250 


-3.00000 

-1.83000 

-3.62250 

3.26266 


0.0 

0.0 

18.00000 

-2.42920 

CROS 

4.52065 

5.19615 

-2.42920 

-4.20750 


-3.39000 

-2.22000 

-4.20750 

3.78955 


# — 

- — 





FINISHED COMPUTATION FOR CROSS SFCT ION 2 
TOTAL NUMBER OF OFF-BODY POINTS GENERATED = 15 


CASE NO. TEST 


PROGRAM J1HP, 


DOUGLAS AIRCRAFT COMPANY 

LONG peach W V1 ?AQH 
THURSDAY, APR 19, 1979 


PAGE 


20 


BODY SECTION NO. 


NO. 


INLET TEST CASE WITH CENTERBODY 1 WAKE 1 PL. SYM. TYPE C 

, 2 TYPE = l TOTAL NO. OF POINTS = 21 

LIFTING SECTION NO. 1 NO. OF SOURCE ELEMENTS 

OF WAKE ELEMENTS 1 TOTAL NO. OF ELEMENTS PER 

TOTAL NO. OF CONTROL POINTS ( INCL. OFF BODY POINTS ) = 


ONE 


NO. 


OF STRIPS 


6 

STRIP 
15 


NO. OF FAR ELEMENTS - 


TOTAL NO. OF ELEMENTS IN THE LIFTING SECTION = 21 

0 NO. OF INTERMEDIATE ELEMENTS = 126 NO. OF 


NEAR ELEMENTS = 


774 


FLOW 

FLUXX 

FLUXY 

FLUXZ 

TCTAL FLUX 

1 

2 

3 

0.9695266*03 
0.5968026*02 
—0. A 29 1536—04 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

8:l?33SSIt84 

-0 . 058307E— 04 

4 

0.4602046*02 

0.0 

0.0 

0.0 

CALCULATION OF CROSS 

SECTION NO. 

1H AS BEEN COMPLETEO 

. 





CASE NO. TEST 


BOOT SECTION NO« 


NO. 


PROGRAM J1HP * DOUGLAS AIRCRAFT COMPANY PAGE 22 

LONG BEACH OIVIS ION 
THURSDAY, APR 19, 1979 

INLET TEST CASE WITH CENT ER30DY 1 WAKE 1 PL. SYM. TYPE C 

3 2 TYPE = 1 TOTAL NO. OF POINTS = 21 

LIFTING SECTION NO. 1 NO. OF SOURCE ELEMENTS 

OF WAKE ELEMENTS l TOTAL NO. OF ELEMENTS PER 

TOTAL NO. OF CONTROL POINTS ( INCL. OFF BODY POINTS ) 
NO. OF ELEMENTS IN THE LIFTING SECTION = 


ONE 


NO. 


FLOW 


NO. OF FAR 
FLUXX 


ELEMENTS = 

FLUXY 


TOTAL 

0 


NO. 
FLUXZ 


OF 


INTERMEDIATE ELEMENTS = 
TCTAL FLUX 


0.920213E+03 0.0 

0. 498369E 4-02 0.0 

0. 3<*3391E-05 0.0 

-0. I04613E+01 0.0 

CALCULATION OF CROSS SECTION NO. 


0 . 0 
0.0 
0.0 
0.0 

2H AS BEEN COMPLETED. 


0 . 184043F4-04 
0. 996 7 37E+0 2 
0. 7867 8 l E— 05 
0.0 


111 


6 

STRIP 7 
= 15 

21 

NO. OF NEAR 


CASE NO. TEST 


PROGRAM J1HP, DOUGLAS AIRCRAFT COMPANY PAGE 

LONG BEACH DIVISION 
THURSDAY, APR 19, 1979 

INLET TEST CASE WITH CENTERBODY 1 WAKE 1 PL. SYM. TYPE C 
FLOW PASSAGE CALCULATION PROGRAM 


23 


ONE 


C.P. 

NO. 


XCP 


CENTERLINE 

YCP 


POINT 


ZCP 


0. 7500GCE+00 0.0 

0. 750000E + 01 0.0 

0*1 50000F 4- 02 0.0 

0.2 10000E 4* 02 0.0 


0.0 

0.0 

0.0 

0.0 


UI 


0 • 1000 00 E 4-0 1 
0. lOOOOOE+Ol 
0. 100000E4-01 
0. 100000E4-01 


NORMAL VECTOR 
UJ 


0.0 

0.0 

0.0 

0.0 


UK 


0.0 

0.0 

0.0 

0.0 


OF STRIPS 


ELEMENTS 


o 

cx> 


NO. i <_.*> i 


f'KUbKAft Jlhr*, 


UuUbL/*^ A 1 K Ltt A t~T 

ILNb b L A Lh 

Phl’JAY, 


1 MLI M.L* wAiL tviTH CtMtKtoL'Y 
I-IiMmL Lin PUI 


WAKi. J P L 


LLftPANY 

u j. VI 5 i Lin 

it* iv7y 


P mU fr 


s y ;-i . iyPl 

uUlVMlwiiULt LbuJl ION 
LuMa i NATION MjLUTaun IiU. 1 


c6 


GNt 


f“ Kc*lL I cU NL 1 I i l.No — 


m:ALc UP h! ( mU\ 
rKU, o I t\£ Art MALn uU- 
1 L l a L I l . ft .P w N A I u 


O • U 


I u 1 A 
fl\LJ 

LImI 

vVlU* 


V 


UtNii I Y 
r\M I i J 
wAL urt-t-L, 




VC. 4 .- 


o • *; j< / tc uc “v^. 

0 • 4 oj U+CJ 


AinClL uF Ya^ 

1 1 J I a L PK^SUKc 
LlAfiC UhFLKAluKc 

LlATlL Ut.NbilY 
I t ft P/b h A — L t V L L KAllu 
MAXIMUM V’rL'uCi l Y 
Lt.UiV oKillUAL VtL. 


t « C' 

^ i +04 

0 .b W Uut+vJ 

0 . ViOt-CL 
J • 1 wCv Jt^Ci 
y • ^'♦V i>^z£ +c *t 
/ oiu + Ci 


HKctilKuAM VlLULITY 
Maul PRcbbkt 
UYKAP.1L PftLbbUKb 


0.1^i600t+U3 

O.«-OYc46i:+04 


l.tNwl I Y KA11U 
PhcL/uL A-lcV*_L KAllU 
pRfchLrfKfc^M LLJ««u LPfctO 
LI ALNA'I iLN LbU^U bPttU 


o.wiwit+co 

O.iuuoOUL + Ojl 
O. l 1 lbY4 t+04 
u.iil^VHt+u^ 


LUlftwL LiAlluiN LbKLlTIONb 


AVLKAww ft A s.n MMtLK 
t v, wrl V A V L K A V" L. VU. 
jfiihMw Kkt^LOrvL 


C • ij/Ov>CC^uC 

u • 1. > A '* Jl I :: + L U 

u. 4 . < i 


li\Lul Mm Lb M. 'JW KAIl u . i 3 l Kwc +0<: 

KL CvlR Pi C • 1 1? b b + ^ 0 

VcLULllY KAliu u • Y <-b -> a 2 t + bu 


LUKK. MAib PLU* RAT b 
KfiiJL UVcK KHuT 
Ai/lKALc VcLlfblTY 


O. JLbbV^Ofc+U*' 
O.VYJbbHb+OO 
O.Ai>33fc>it+OJ 


otftru 1 c 


n m w b ft Uj W 


0 . / V <w Itc + O 
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LONG BEACH OIVISION 
CASE NO. TEST THURSDAY, APR 19, 1979 


TNLET TEST 

CASE WITH CENTERBODY 1 WAKE 1 PL 

. SYM. TYPE 

C ONE 


FINAL OUTPUT 

COMPRESSIBLE SOLUTION 



NON-SYMMETRY PRINT OUTPUT 



SOLUTION 

FOR THE ON-BODY 

ELEMENTS 


N M 

xc 

YC 

ZC 

VX 

VY 

vz 

VM AG 
MACH 
P / PT 

VBAR 
EPSLON 
AREAtFP ) 

1 1 

7.999998 

-0.346410 

0.600000 

255.825241 

-44.310364 

43.679840 

263.282959 

0.237253 

0.961577 

266.170654 
C. 970045 
112.236 832 

2 

10.636363 

-0.732185 

1.268181 

33 7.686768 
-43.867065 
42.717133 

343. 193115 
0.310487 
0.935346 

321.380371 
C. 955167 
92.955 750 

3 

18.000000 

-0.909327 

1.574998 

342.087402 

-0.000198 

2.538933 

342.096924 

0.309476 

0.935750 

364.805664 

0.940721 

81.890594 

************ 


************ 


************ 

2 1 

7.999998 
-0.692819 
0. JOOOOO 

256.037109 

-88.693588 

-0.000158 

270.964355 

0.244257 

0.959337 

266.170654 
C. 970045 
112.236832 

2 

10.636362 

-1.464368 

-0.0000.00 

337.948975 

-87.801758 

0.000023 

349.168701 

0.316000 

0.933132 

321.380371 

0.955167 

92.955750 

3 

18.000000 

-1.818653 

0.0 

342.091553 

0.000015 

0.000036 

342.091797 

0.309472 

0.935750 

364.873535 
C. 940 69 6 
81.875381 

************ 


************ 


************ 

3 1 

7.999998 

-0.346410 

-0.600000 

255.825653 

-44.309662 

-43.679916 

263.283447 

0.237254 

0.961577 

266.170654 

0.970045 

112.236832 

2 

10.636363 

-0.732185 

-1.263180 

337.686768 

-43.866455 

-42.717224 

343. 193115 
0.310487 
0.935346 

321.380371 

0.955167 

92.955750 

3 

18.000000 

-0.909326 

-1.574998 

342.090820 

0.000070 

-2.538962 

342.100342 
C. 309480 
0.935747 

364.856445 

0.940703 

81.879181 

************ 


************ 


************ 

4 1 

18.000000 

-2.598071 

4.500000 

337.698975 

-0.000143 

2.423266 

337.707764 

0.305431 

0.937352 

365.448975 

0.940487 

81.746475 

2 

7.500000 

-2.857880 

4.949999 

305.916016 

17.661087 

-14.202209 

306.754395 

0.276985 

0.948100 

266.170654 
C. 970045 
112.236816 

3 

1.500000 

-3.767204 

6.524998 

102.800751 

44.513550 

-36.211075 

117. 731537 
0. 105617 
0.992232 

143.946915 
C. 9,91 560 
207.536011 

4 

1.500000 

-5.066241 

8.774994 

121.370392 
-52.55493 2 
97.876053 

164.537415 

0.147764 

0.984868 

56.035858 
G. 998737 
533.125732 

5 

7.500000 

-5.975566 

10.349991 

165.750153 

-9.569506 

20.959061 

167. 344025 
0. 150296 
0.984348 

16.703252 

0.999888 

1787.98730 

6 

18.000000 

-6.235375 

10.799988 

135.253372 

-0.000207 

4.279458 

135.321167 
0. 121441 
0.989747 

5.673552 

0.999987 

5265.51562 

************ 


*** ********* 


****** ****** 
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PROGRAM JIHP* DOUGLAS AIRCRAFT COMPANY PAGE 29 

LONG BEACH OIVISION 
CASE NO. TEST THURSDAY, APR 19, 1979 

INLET TEST CASE WITH CENTERBQOY 1 WAKE l PL. SYM. TYPE C CNE 

FINAL OUTPUT COMPRESSIBLE SOLUTION 

NON-SYMMETRY PRINT OUTPUT 




SOLUTION 

FOR THE ON-BODY 

ELEMENTS 


N 

M 

W 

zc 

vx 

VY 

VZ 

VM AG 
MACH 
P / PT 

E^SLON 
AREAIFP I 

5 

1 

18.000000 

-5.196150 

0.0 

33 7.69848 6 
-0.000251 
0.000008 

337.698730 

0.305423 

0.937355 

365.448975 
0.940 48 7 
81. 746475 


2 

7.500000 

-5.715766 

0.0 

305.956543 

35.328491 

0.000058 

307.989746 

0.270117 

0.947689 

266.170654 
C. 970045 
112.236816 


3 

1.500000 

-7.534418 

0.0 

102.910522 

89.123047 

0.000048 

136. 137894 
0. 1221 76 
0.989622 

143.946915 
C. 991 560 
207.536011 


4 

1.500000 

-10.132492 

0.0 

121.359207 

-105.100464 

0.000050 

160.543488 
0. 144162 
0.985589 

56.035858 

0.998737 

533.125732 


5 

7.500000 

-11.951141 

0.0 

165.749908 

-19.139160 

0.000033 

16 6. 851410 
0. 149852 
0.984440 

16.708252 
0.999888 
L787. 98730 


6 

18.000000 

-12.470764 

0.0 

135.253296 

-0.000963 

0.000012 

135.253403 

0.121379 

0.989758 

5.673552 

0.999987 

5265.51562 

************ 


************ 


************ 

6 

1 

18.000000 

-2.598075 

-4.500000 

337.698730 

-0.000025 

-2.423203 

337.707764 

0.305431 

0.937352 

365.448975 

0.9404B7 

81.746475 


2 

7.500000 

-2.857882 

-4.949999 

305.916016 

17.661621 

14.202299 

306.754395 

0.276985 

0.948100 

266.170654 
C.970 045 
112.236816 


3 

1.500000 

-3.767209 

-6.524998 

102.800186 

44.513626 

36.211151 

117.731094 
0. 105617 
0.992232 

143.946915 

0.991560 

207.536011 


4 

1.500000 

-5.066248 

-8.774994 

121.370255 

-52.554836 

-97.876063 

164.537323 

0.147764 

0.984868 

56.035858 

0.990737 

533.125732 


5 

-10.349991 

16 5.750320 
-9.569611 
-20.959091 

167.344193 

0.150296 

0.984348 

16.700252 

0.999888 

1787.98730 


6 

18.000000 

-6.235385 

-10.799988 

135.253418 

-0.000322 

-4.279426 

135. 321220 
0. 121441 
0.989747 

5.673552 

0.999987 

5265.51562 

************ 


************ 


************ 
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PROGRAM JIHP, DOUGLAS AIRCRAFT COMPANY PAGE 31 

LONG BEACH OIVISION 
CASE NO. TEST THURSDAY, APR 19, 1979 

INLET TEST CASE WITH CENTERBODY l WAKE l PL. SYM. TYPE C ONE 




FINAL 

OUT PUT 

COMPRESS IDLE SOLUTION 




SYMMETRY 

PRINT 

OUTPUT — REFLECTED ELEMENTS 




SOLUTION 

FOR THE ON- BODY ELEMENTS 


N 

M 

XC 


VX VMAG 

VBAR 



YC 


VY MACH 

EPSLON 



ZC 


VZ P / PT 

AR EA ( FP ) 


1R l 

7.99999ft 

255.805099 

262.551758 

266.170654 


0.346410 

44.306885 

0.236587 

0.970045 


0.600000 

39.166000 

0.961787 

112.236832 

2 

10.636363 

337.664795 

342.695068 

321.380371 


0.732185 

43.864227 

0.310028 

0.955167 


1 . 268 L81 

38.705170 

0.935529 

92.955750 

3 

18.000000 

342.087402 

342.096924 

364.805664 


0.909327 

0.000198 

0.309476 

C. 940 721 


1.574998 

-2.538939 

0.935750 

81.890594 

************ 


************ 


************ 

2R l 

7.999998 

256.037109 

270.964355 

266.170654 


0.692819 

88.693588 

0.244257 

C. 970045 


0.000000 

-0.000133 

0.959337 

112.236832 

2 

10.636362 

337.948975 

349.168701 

321.380371 


1.464368 

87.801758 

0.316000 

0.955167 


- 0.000000 

0.000021 

0.933132 

92.955750 

3 

18.000000 

342.091553 

342.091797 

364.873535 


1.818653 

-0.000015 

0.309472 

C.940696 


0.0 

0.000042 

0.935750 

81.875381 

************ 


************ 


************ 

3R 1 

7.999998 

255.805511 

262.552002 

266.170654 


0.346410 

44.306168 

0.236587 

C. 9 70045 


-0.600000 

-39. 166046 

0.961787 

112.236832 

2 

10.636363 

337.664795 

342.694824 

321.380371 


0.732185 

43.863617 

0.310028 

0.955167 


-1.268180 

-38.705215 

0.935529 

92.955750 

3 

18.000000 

342.090820 

342.100342 

364.856445 


0.909326 

-0.000070 

0.309480 

C. 940 703 


-1.574998 

2.538970 

0.935747 

81.879181 

************ 


************ 


************ 

4R 1 

18.000000 

33 7.698975 

337.707764 

365.448975 


2.598071 

0.000143 

0.305431 

0.940487 


4.500000 

-2.423064 

0.937352 

81.746475 

2 

W 

305.924072 

307.004639 

266.1 70654 


-17.661560 

0.277214 

0.970045 


4.949999 

-18.716019 

0.948016 

112.23681 6 

3 

1 . 500000 

102.799973 

117.609390 

143.946915 


3.767204 

-44.513214 

0. 105507 

0.991560 


6.524998 

-35.814575 

0.992248 

207.536011 

4 

1.500000 

121.370621 

164.622757 

56.035858 


5.066241 

52.555023 

0. 147841 

0.998 737 


8.774994 

98.019104 

0.984852 

533.125732 

5 

7.500000 

165.749207 

166.508484 

16.708252 


5.975566 

9.569451 

0 . 149542 

0.999888 


10.349991 

12.674822 

0.984505 

1787.98730 

6 

18.000000 

135.253372 

135.321167 

5.673552 


6.235375 

0.000207 

0.121441 

0.990987 


10.799988 

-4.279350 

0.989747 

5265.51562 


************ ****** 4 ***** ************ 
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PROGRAM J1HP, DOUGLAS AIRCRAFT COMPANY PAGE 33 

LONG BEACH DIVISION 
CASE NO. TEST THURSDAY, APR 19, 1979 

INLET TEST CASE WITH CENTERBODY 1 WAKE l PL. SYM. TYPE C ONE 



FINAL OUTPUT 

COMPRESSIBLE SOLUTION 



SYMMETRY PRINT 

OUTPUT — REFLECTED ELEMENTS 



SOLUTION 

FOR THE ON- BODY 

ELEMENTS 


M 

xc 

VX 

VM AG 

VBAR 


YC 

VY 

MACH 

EPSLON 


ZC 

VZ 

P / PT 

AREAIFP ) 

1 

18.000000 

33 7.698486 

337.698730 

365.448975 


5.196150 

0.000251 

0.305423 

0.940487 


0.0 

0.000016 

0.937355 

81.746475 

2 

7.500000 

305.956543 

307.989746 

266.170654 


5.715766 

-35.328491 

0.273117 

C. 970045 


0.0 

0.000064 

0.9476B9 

112.236816 

3 

1.500000 

102.910522 

136. 137894 

143.946915 


7.534418 

-89.123047 

0.122176 

C. 991560 


0.0 

0.000050 

0.989622 

207.536011 

4 

1. 500000 

121.359207 

160.54348 8 

56.035858 


lg. ^32492 

105.100464 

0.000053 

0. 14416? 
0.985589 

C * 9 9 8 7 3 7 
533.125732 

5 

7.500000 

165.749903 

166.851410 

16.708252 


11.951141 

19.139160 

0.149852 

0.999888 


0.0 

0.000030 

0.984440 

1787.98730 

6 

18.000000 

135.253296 

135.253403 

5.673552 


12.470764 

0.000963 

0. 121379 

0.999987 


0.0 

0.000008 

0.989758 

5265.51562 


************ ************ ************ 


6R 1 

18.000000 

2.598075 

-4.500000 

337.698730 

0.000025 

2.423139 

337.707764 

0.305431 

0.937352 

365.448975 

0.940487 

81.746475 

2 

7.500000 

2.857882 

-4.949999 

305.924072 

-17.662094 

18.716110 

307.004639 

0.277214 

0.948016 

266.170654 

0.970045 

112.236816 

3 

1.500000 

3*767209 

-6.524998 

102.799408 

-44.513290 

35.814651 

117.608932 

0.105507 

0.992248 

143.946915 

0.991560 

207.536011 

4 

1.500000 

5.066248 

-8.774994 

121.370483 

52.554977 

-98.019104 

164.622650 

0.147841 

0.984852 

56.035850 

0.998737 

533.125732 

5 

7.500000 

5.975574 

-10.349991 

165.749374 

9.569556 

-12.674834 

166.508652 
0. 149542 
0.984505 

16.703252 

0.999888 

1787.98730 

6 

18.000000 

6.235385 

-10.799988 

135.253418 

0.000322 

4.279392 

135.321228 

0.121441 

0.989747 

5.673552 

0.999987 

5265.51562 

************ 


************ 


************ 
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CASE 

PROGRAM 

NO. TEST 
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43 



INLET 

TEST CASE WITH CENTERBODY 1 WAKE 

1 PL. SYM. TYPE C 

ONE 





TABLE OF RUNT IME 

BREAKDOWN 



SUBROUTINES 

CPU TIME 
AT CALL 
STATEMENT 

CPU TIME 
AFTER RETURN 
TO CALL PROG 

CPU TIME 
USED 

I/O TIME 
AT CALL 
STATEMENT 

I/O TIME 
AFTER RETURN 
TO CALL PROG 

I/O TIME 
USED 

MAIN 

0.000 

16.590 

16.590 

0.0 

13.244 

13.244 

REPAN 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

INPUT 

0.035 

0. 632 

0.597 

0.068 

1.3 86 

1.318 

V FORM 

0.632 

2.194 

1.562 

1.386 

2.286 

0.900 

PSWISE 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

AFORM 

2. 104 

2.564 

0.370 

2.286 

2.566 

0.280 

COL SOL 

2.564 

3.233 

0.668 

2.566 

3.578 

1.012 

FSLUTN 

4.100 

4.824 

0.724 

3.376 

4.066 

0.690 

CROSEC 

4. 395 

7.853 

2.958 

4.186 

7.488 

3.302 

FLPASG 

7.886 

8.501 

0.615 

7.408 

8.804 

1.316 

COM PRO 

0.0 

0. 0 

0.0 

0.0 

0.0 

0. 0 

ICMPRS 

8.501 

16.590 

8.088 

8.804 

13.244 

4.440 


( UNIT IN SECONDS ) 

*** END OF INPUT DATA. NORMAL PROGRAM TERMINATION. 
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